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Abstract 

The  objective  of  this  research  was  to  determine  the  feasibility  of  using  a  nuclear 
reaction  heat  source,  such  as  the  electromagnetically  triggered  decay  of  an  isomer, 
in  a  solid-state  heat  exchanger  to  power  an  off-the-shelf  gas  turbine  engine.  Two 
primary  performance  measures  examined  were  the  total  pressure  decrement  across 
the  heat  exchanger  and  the  total  temperature  capability  leaving  the  heat  exchanger. 

The  analysis  included  the  use  of  a  commercial  software  package,  ANSYS®  5.6.1, 
running  on  a  700  MHz  Pentium®  III  PC.  This  package  includes  the  FLOTRAN® 
computational  fluid  dynamics  program,  a  finite  element  program  based  on  unstruc¬ 
tured  meshes,  with  multiple  discretization  schemes,  turbulence  models,  and  advec- 
tion  options.  Boundary  conditions  on  velocity,  pressure,  temperature,  heat  flux,  and 
heat  generation  are  available  and  were  used  in  this  research. 

Three  basic  geometries  of  heat  exchanger  were  explored  in  this  research:  con¬ 
centric  annular  tubes,  radial  trapezoidal  fins,  and  a  dual,  concentric  annulus  of  rect¬ 
angular  fins.  These  were  selected  due  to  the  simplicity  of  geometry  and  potential 
ease  of  manufacture.  In  addition,  because  the  flow  through  all  of  these  geometries 
could  be  reasonably  approximated  by  a  series  of  two-dimensional  flow  fields,  run 
times  were  on  the  order  of  1  day,  a  significant  reduction  from  3-D  flow  calculations. 

All  three  configurations  produced  sufficient  heat  transfer.  Pressure  ratios 
across  the  heat  exchangers  varied  in  the  range  from  94.5%  to  97.5%.  Turbine  inlet 
temperatures  varied  from  986  K  to  1150  K  (1775°R  to  2070°R).  In  the  J-57  engine, 
these  conditions  will  produce  a  static,  sea-level  thrust  of  approximately  37,000  N 
(8,300  lb.)  to  47,000  N  (10,600  lb.),  compared  to  46,000  N  (10,300  lb.)  for  the 
conventional  engine. 
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ANALYSIS  OF  THE  APPLICATION  OF  A  TRIGGERED  ISOMER 
HEAT  EXCHANGER  AS  A  REPLACEMENT  FOR  THE 
COMBUSTION  CHAMBER  IN  AN  OFF-THE-SHELF  TURBOJET 


I.  Introduction 

1.1  Motivation 

Since  1942,  with  the  creation  of  the  first  fission  chain  reactor,  scientists  and 
engineers  have  contemplated  the  use  of  atomic  energy  to  power  aerospace  vehicles. 
These  proposals  began  with  the  idea  of  flowing  light  fluids  through  a  nuclear  reactor 
to  power  a  rocket  motor.  Starting  in  1946,  and  continuing  through  the  1950s  and 
1960s,  the  United  States  Air  Force  and  the  Atomic  Energy  Commission  worked  on 
projects  including  Nuclear  Propulsion  of  Aircraft  (NEPA)  and  later.  Aircraft  Nuclear 
Propulsion  (ANP),  which  explored  the  use  of  nuclear  reactors  to  power  jet  aircraft. 
Ultimately,  these  projects  were  cancelled  due  to  a  technical  inability  to  produce  a 
safe,  controllable  reactor  that  met  the  power  requirements  of  flight  at  a  sufficiently 
low  weight,  as  well  as  safety  concerns  for  the  crew  and  civilians  who  might  be  affected 
(1:1-6) 

Recent  developments  in  controlling  the  decay  of  several  long  half-life  isomers 
of  hafnium,  lutetium,  and  tantalum,  may  have  opened  another  possibility.  Some  of 
these  isomers  have  extremely  high  internal  energy  states,  and  an  atomic  structure 
which  inhibits  rapid  decay  to  lower  energy  states.  By  triggering  rapid  decay  of  these 
isomers  through  X-ray  bombardment,  a  controllable  release  of  the  internal  atomic 
energy  available  may  be  possible.  The  specific  isomer  of  interest  for  this  study  was  the 
4-quasiparticle  isomer  of  hafnium,  which  has  a  half-life  of  31  years  and  an  excitation 
energy  of  2.446  mega-electron  volts  (MeV).  This  has  the  potential  to  release  up  to 
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1.3  gigajoules  per  gram,  mostly  in  the  form  of  heat  and  some  gamma  rays.  This 
could  represent  an  extremely  compact,  controllable  heat  source  with  the  potential  to 
power  an  aircraft  engine  at  a  fuel  consumption  rate  on  the  order  of  ten  pounds  per 
day  (2;695;3:1).  This  power  source  may  represent  a  practical  method  of  achieving 
extraordinary  endurance,  free  of  many  of  the  safety  and  environmental  issues  which 
plagued  fission-powered  flight  programs  in  the  past. 

Practical  applications  of  this  type  of  power  source  could  be  widespread,  de¬ 
pending  on  the  expense  and  shielding  requirements,  but  long-endurance,  unmanned 
aircraft  are  one  of  the  first  obvious  applications.  If  the  expense  can  be  reduced  suffi¬ 
ciently,  this  power  source  could  be  practical  for  any  application  where  long  endurance 
is  a  desirable  characteristic.  In  any  application,  the  additional  weight  of  radiation 
shielding  over  any  sensitive  components  of  the  aircraft  (including  personnel)  would 
have  to  be  taken  into  account. 

1.2  Problem  Statement 

The  goal  of  this  study  was  to  determine  if  a  solid-state  heat  exchanger  con¬ 
figuration  could  be  devised  that  could  feasibly  replace  the  combustion  chamber  in 
an  off-the-shelf  gas  turbine  engine.  Solid-state  is  taken  to  mean  a  heat  exchanger 
that  does  not  rely  on  a  fluid  for  the  heat  source.  Thus,  the  heat  exchanger  will 
generate  its  own  heat,  rather  than  transferring  heat  from  a  hot  fluid  to  a  cool  fluid. 
For  the  heat  exchanger  to  be  feasible,  it  must  meet  several  requirements.  First,  the 
flow  through  the  heat  exchanger  must  reach  stagnation  temperatures  sufficient  to 
produce  thrust.  Second,  the  stagnation  pressure  losses  through  the  heat  exchanger 
must  be  low  enough  to  allow  operation  of  the  engine.  In  addition,  the  added  weight 
of  the  heat  exchanger  and  the  necessary  radiation  shielding  around  sensitive  areas 
must  not  be  so  high  as  to  require  extensive  re-design  of  the  airframe  to  accommodate 
the  engine.  Last,  the  heat  exchanger  should  fit  within  the  existing  engine  envelope 
without  requiring  extensive  mechanical  changes  to  the  engine  systems  such  as  the 
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drive  spools  for  the  compressor  stages.  To  some  extent,  all  of  these  criteria  are  appli¬ 
cation  driven.  The  temperature  gain  and  pressure  loss  must  provide  adequate  thrust, 
but  exactly  what  level  of  thrust  is  adequate  is  driven  by  the  intended  application, 
including  airframe  design  and  mission.  Therefore,  for  now,  it  was  decided  to  address 
the  adequacy  of  the  thrust  level  only  as  a  percentage  of  the  original  design  thrust  of 
the  chosen  gas  turbine  engine  (4).  It  was  decided  to  address  the  weight  issues  only 
in  general  terms,  as  they  will  be  driven  almost  entirely  by  the  precise  application. 
For  example,  how  much  radiation  shielding  is  required  will  depend  heavily  on  the 
airframe  geometry  and  the  amount  of  radiation  that  the  most  sensitive  components 
can  tolerate,  while  the  additional  weight  of  the  heat  exchanger  will  depend  on  factors 
including  (but  not  limited  too):  the  air  loads  on  the  fins,  the  design  mission  of  the 
aircraft  and  engine,  the  turbine  inlet  temperature  required,  and  the  mass  flow  of  air 
through  the  heat  exchanger. 

1.3  Summary  of  Current  Knowledge 

Current  knowledge  in  this  type  of  heat  exchanger  can  be  divided  into  two 
basic  areas:  work  done  in  the  past  with  aircraft  and  rockets  to  be  powered  by  fission 
reactions,  and  work  done  with  more  conventional  heat  exchangers  and  augmenting 
the  heat  transfer. 

Most  work  with  nuclear  powered  aircraft  was  accomplished  in  the  1950s  and 
1960s,  and  concentrated  on  two  concepts.  Both  used  fission  reactors,  but  differed  in 
the  mechanism  used  to  heat  the  air.  One  was  a  direct  cycle,  in  which  the  air  was 
passed  directly  through  the  reactor  core,  located  in  the  engine.  The  second  was  an 
indirect  cycle,  in  which  the  reactor  heated  a  liquid  metal,  which  was  circulated  to  a 
liquid-to-air  heat  exchanger  in  the  engine  (1:3). 

Of  these  two,  the  first  was  most  directly  applicable  to  the  current  work,  as  it  in¬ 
volves  solid-state  heat  exchangers,  without  tubes  to  significantly  disturb  the  airflow. 
Various  simple  configurations  were  explored,  in  terms  of  advantages  and  disadvantage 
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in  the  fission  reactor  application.  These  included  flat  plates,  solid  rods,  and  tubes. 
These  are  shown  schematically  in  Figure  1.1.  Flat  plates  hold  the  advantages  of  ease 
of  manufacture,  easy  structural  support  and  a  near  immunity  to  blocked  passages. 
Debris  and  uneven  fuel  loading  will  not  drastically  affect  the  flow  in  a  channel,  as 
the  flow  will  simply  go  around  the  area.  Solid  rods  have  similar  advantages  to  plates, 
except  that  structural  support  is  significantly  more  complicated,  and  non-uniform 
distribution  of  heat-generating  material  can  lead  to  significant  temperature  gradi¬ 
ents.  Drilling  holes  in  a  solid  block,  or  stacking  interlocking  cylinders  would  lead  to 
the  tube  configuration.  This  is  the  most  structurally  stable  configuration,  however, 
it  is  also  the  most  subject  to  blocking  and  failure  due  to  local  hot  spots  (1:152;5;6). 
Except  for  the  difference  in  materials  (uranium  instead  of  hafnium)  and  the  concur¬ 
rent  difference  in  energy  content,  this  is  a  close  analog  of  the  problem  at  hand.  In 
the  interest  of  simplicity  of  manufacturing  and  support,  all  of  the  heat  exchanger 
configurations  used  in  this  research  were  variations  on  the  flat  plate  configuration. 


Figure  1.1  Basic  Configurations  of  Heat  Exchanger 
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For  the  purposes  of  this  study,  the  second  method  of  transferring  heat  to 
the  air  from  the  nuclear  reactor  overlaps  with  the  research  on  conventional  heat 
exchangers.  In  both  cases,  heat  is  transferred  from  a  liquid  traveling  in  tubes  through 
fins  attached  to  the  tubes,  the  difference  being  in  the  nature  of  the  liquid,  with  the 
reactor  using  liquid  metals,  such  as  sodium.  The  aim  with  most  heat  exchangers  of 
this  type  is  to  cool  the  liquid,  rather  than  heat  the  air,  leading  to  less  attention  being 
paid  to  the  efficiency  of  the  air  side  of  the  heat  exchanger,  at  least  with  respect  to 
pressure  losses.  Much  of  the  effort  in  increasing  heat  transfer  is  in  various  methods 
of  augmenting  the  heat  transfer,  usually  by  inducing  turbulence.  Devices  such  as 
ribbed  fins,  interrupted  fins,  and  roughened  surfaces  are  used  to  create  turbulence 
and  higher  heat  transfer  to  the  air  at  the  cost  of  reduced  total  pressure  across  the 
heat  exchanger.  In  the  liquid  side  of  the  liquid-air  heat  exchanger,  twisted  ribbons 
and  similar  devices  can  be  used  to  increase  the  heat  transfer  from  the  liquid  (also  by 
inducing  additional  turbulence).  In  some  heat  exchangers,  multiple  passes  for  one 
or  both  fluids  are  used  to  increase  the  heat  transfer  from  a  given  volume  of  fluid 
(7;8;9:417-477;10). 

There  are  also  active  methods  for  augmenting  heat  transfer,  including  vibration 
of  the  surfaces  of  fluids,  electrostatic  fields,  and  suction  or  injection  jets  onto  the 
surfaces,  but  these  would  add  significant  weight  to  an  already  hefty  gas  turbine 
engine  and  hafnium  heat  exchanger  combination  (7). 

1.4  Present  Work 

The  present  work  used  the  geometry  and  combustion  chamber  inlet  condi¬ 
tions  of  the  J-57  turbojet,  used  in  the  Boeing  707,  KC-135  and  some  models  of 
B-52.  This  engine  was  chosen  largely  for  its  simple  geometry  and  the  availability  of 
pressure,  temperature,  and  geometric  data.  This  engine  is  shown  schematically  in 
Figure  1.2.  As  shown  in  the  figure,  heat  exchangers  (manufactured  from  hafnium 
isomer)  would  be  installed  in  place  of  the  combustion  chamber.  All  of  the  heated 
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surfaces  of  these  heat  exchangers  were  assumed  to  be  made  of  the  hafnium  isomer. 
Therefore,  the  heat  exchanger  itself  was  also  the  heat  source.  Simulations  using 
ANSYS® /FLOTRAN^^  (version  5.6.1)  (11),  on  a  Pentium®  III  platform,  were  run 
to  evaluate  the  effectiveness  of  several  heat  exchanger  configurations  in  transferring 
sufficient  heat  to  the  airflow  with  low  pressure  loss.  These  heat  exchanger  configura¬ 
tions  were  designed  under  the  assumption  that  the  isomer  heat  source  would  make 
up  the  heated  surfaces  of  the  heat  exchanger.  Since  this  was  only  a  feasibility  study, 
very  little  has  been  done  to  optimize  the  configurations,  or  attempt  to  exceed  the 
performance  of  the  original  engine. 


Combustion  Chamber/ 
Heat  Exchanger 


High-Low  Pressure 


Nozzle 


Air  Flow 


Figure  1.2  Schematic  View  of  J-57  Turbojet 
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1.5  Outline 


This  section  contains  an  overview  of  the  chapters  to  follow.  Chapter  2  con¬ 
tains  a  discussion  of  the  theory  and  equations  which  pertain  to  the  current  work. 
Discussion  of  the  methods,  hardware,  and  software,  including  inputs  to  the  software 
(geometries,  grids,  and  boundary  conditions)  used  in  this  work  follows  in  Chapter  3. 
The  results  of  computation  are  presented  in  Chapter  4,  including  some  estimates  of 
performance  relative  to  the  baseline  (petroleum  powered)  engine.  Conclusions  and 
recommendations  for  future  research  in  this  area  are  contained  in  Chapter  5. 
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11.  Background  and  Theory 

This  chapter  contains  a  discussion  of  the  Navier-Stokes  equations  for  fluid 
flow,  heat  transfer,  and  the  turbulence  model  used  in  the  computation,  as  well  as 
background  information  on  the  process  of  triggered  hafnium  decay  and  estimating 
the  radiation  shielding  requirements. 


2.1  Navier-Stokes  Equations 

The  program  used  in  this  research,  ANSYS®  (11),  implements  the  Navier- 
Stokes  equations  through  a  finite-element  method  of  approximating  differentials.  To 
understand  the  operation  of  this  program,  it  is  necessary  to  understand  the  basic 
equations  on  which  it  is  built.  Therefore,  this  section  includes  a  basic  presentation 
and  explanation  of  the  Navier-Stokes  equations. 


2.1.1  Substantial  Derivative  Operator.  Before  beginning  a  presentation  of 
the  Navier-Stokes  equations,  it  is  imperative  that  the  reader  understand  the  ^  op¬ 
erator  for  the  substantial  derivative.  This  allows  a  much  more  compact  presentation 
of  the  equations  than  would  otherwise  be  possible.  The  ^  operator  can  be  defined 


as  (12:61-73;13): 


£0 

Dt 


d()  d{)  dO  d{) 

dt  dx  dy  dz 


(2.1) 


and  will  be  used  repeatedly  in  this  presentation. 


2.1.2  Continuity  Equation.  If  the  conservation  of  mass  law  is  applied  to  a 
fluid  passing  through  an  infinitesimal,  fixed  control  volume,  the  continuity  equation 
can  be  derived: 

^  +  V-{pf)  =  0  (2.2) 

where  p  is  the  fluid  density  and  is  the  vector  fluid  velocity.  The  first  term 
represents  the  rate  of  change  of  density  in  the  control  volume,  and  the  second  term 
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represents  the  rate  of  mass  flow  out  of  the  control  volume  per  unit  volume  (12:61- 
73;14:3-6;15:250).  This  can  be  further  reduced  in  some  cases  to  eliminate  changes 
in  density,  but  not  in  this  case.  Since  the  problem  was  steady  state,  it  would  seem 
that  the  time  derivative  term  could  have  been  dropped,  but  this  was  not  the  case. 
The  nature  of  the  equation  set  is  such  that  the  problem  is  more  easily  addressed  in  a 
transient  sense,  where  finite  time  derivatives  are  allowed.  In  a  steady-state  problem, 
this  equation  set  is  elliptic,  where  conditions  at  any  location  affect  conditions  at  all 
locations.  This  is  computationally  more  difficult.  By  making  the  problem  transient, 
the  equation  set  is  hyperbolic,  and  computationally  easier  (13). 


2.1.3  Momentum  Equation.  By  applying  Newton’s  second  law  to  a  similar 
control  volume,  and  simplifying  the  result  using  the  continuity  equation  (Equation 
2.2),  the  following  momentum  equations  can  be  derived  (in  the  x,  y,  and  z  directions, 
respectively): 
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(2.3) 


where  fx,  fy,  and  are  the  body  forces  in  the  x,  y,  and  z  directions,  p  is  the  absolute 
viscosity,  and  u,  v,  and  w  are  the  components  of  velocity  in  the  x,  y,  and  z  directions. 
This  formulation  ignores  bulk  viscosity  and  assumes  that  we  have  a  Newtonian  fluid. 
For  the  purposes  of  the  problem  covered  here,  body  forces  were  neglected,  as  the 
problem  is  largely  at  a  single  altitude,  and  the  only  body  force  expected  would  have 
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been  gravity.  Again,  the  problem  was  treated  as  transient,  and  fi  was  left  inside  the 
differential,  since  large  temperature  changes  are  expected,  and  viscosity  is  a  function 
of  temperature  (12;61-73;14:3-6;15:253) 


2.1.4  Energy  Equation.  By  applying  the  First  Law  of  Thermodynamics  to 
the  fluid  passing  through  the  same  control  volume,  and  neglecting  shaft  work,  the 
following  energy  equation  was  found: 


where  h  is  enthalpy,  P  is  pressure,  Q  is  radiative  heat  transfer,  is  conductive  heat 
transfer,  and  the  terms  in  brackets  make  up  the  viscous  dissipation  element  (12:61- 
73;14:3-6;15:257).  The  challenge  at  this  point,  for  this  research  was  to  define  the 
Q  is  neglected  and,  once  again,  the  transient  terms  are  kept  for  computational 
purposes. 

2.1.5  Perfect  Gas  Relationship.  To  make  a  complete  equation  set  of  the 
Navier-Stokes  relations  (Equations  2.2,  2.3,  and  2.4)  an  additional  equation  is  nec¬ 
essary  (assuming  closed  relationships  are  known  for  properties  such  as  viscosity, 
conductivity  and  so  forth).  This  is  normally  the  perfect  gas  relationship,  seen  below: 

P  =  pRT  (2.5) 

where  P  is  pressure,  p  is  density,  T  is  temperature,  and  R  is  the  specific  gas  constant 
for  the  fluid  involved.  R  can  be  related  to  the  universal  gas  constant  by  the  following 
formula: 

R=h  (2.6) 

M 
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where  Ru  is  the  universal  gas  constant,  and  M  is  the  molar  mass  of  the  gas  in 
question.  This  gives  a  complete  set  of  equations,  relating  density  to  local  pressure 
and  temperature  at  any  location  (12:40;15:258). 

2.2  Turbulence  Model 

In  the  process  of  constructing  a  computational  fluid  dynamic  model,  a  choice 
must  be  made  about  how  to  handle  turbulence.  If  suflElcient  computer  power  were 
available,  an  incredibly  fine  mesh  might  be  established,  and  no  modeling  would 
be  necessary.  This  would  require  a  mesh  of  hundreds  of  thousands  or  millions  of 
elements  for  even  the  smallest  problems,  so  this  is  impractical.  One  of  the  best 
accepted  alternatives  is  a  2  equation  model  loosely  based  on  the  Bousinesq  mixing 
length  hypothesis  (16;17;18:42).  This  hypothesis  provides  a  model  for  turbulent 
viscosity,  Ht,  where: 

(2-7) 

oy 

The  problem  is  then  to  quantify  the  mixing  length,  /.  Shih,  Zhu,  and  Lumley  pro¬ 
vide  a  model  in  their  paper  (19)  which  uses  2  new  partial  diflferential  equations  for 
turbulent  kinetic  energy  and  turbulent  dissipation  to  supplement  the  Navier-Stokes 
equations.  In  this  case  fj,t  is  defined  through  the  following  relation: 

IJ^t  =  (2-8) 

where  k  and  e  are  calculated  from  partial  differential  equations  of  the  form: 

^=\(p+^)k]-WhUij-e  (2.9) 

Dt  L  \  J  i 

and 
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where  Cei,  0^2,  c^k  and  are  all  constants,  specified  in  the  model.  This  system  allows 
the  laminar  stresses  (due  to  the  dynamic  viscosity  and  the  mean  velocity  gradients) 
and  the  turbulent  stresses  (modeled  using  and  the  mean  velocity  gradients)  to 
be  combined  into  a  single  stress,  which  can  be  calculated  for  each  element  in  the 
flowfield.  Using  this  system,  the  overall  behavior  of  the  flowfield  can  be  modeled, 
including  some  of  the  turbulent  effects,  without  excessive  computational  require¬ 
ments.  Since  this  did  not  exactly  represent  all  of  the  possible  loss  mechanisms,  it 
cannot  capture  all  of  the  losses,  but  it  should  have  captured  the  large  majority  of 
the  losses  and  allowed  more  rapid  determination  of  feasibility  for  the  heat  exchanger 
configurations  under  consideration  (19). 

2.S  Heat  Transfer 

2.3.1  Heat  Flow  in  the  Heat  Exchanger.  Since  the  heat  exchanger  was 
assumed  to  be  manufactured  from  the  hafnium  isomer,  the  heat  is  generated  within 
the  solid  isomer  material  of  the  heat  exchanger.  This  is  assumed  to  occur  at  a 
constant  volumetric  rate  (^).  Heat  would  then  be  convected  off  of  the  surface,  or 
conducted  along  the  surface.  The  heat  flow  in  a  control  volume  is  shown  in  Figure 
2.1.  In  steady  state  operation,  the  total  of  the  heat  fluxes  on  a  control  volume  must 
be  zero  (i.e.  qgenerated  +  qc^d-ir.  -  qcond.out  -  Qconv.  =  0).  Since  the  q  for  conduction  is 
given  by: 


qcond.  =  k—t  (2.11) 

and  ^generated  1®  given  by. 

^generated  =  {HGR)dV  (2.12) 

where  k  is  thermal  conductivity,  t  is  thickness,  ^  is  the  temperature  gradient,  HGR 
is  the  volumetric  heat  generation  rate,  and  dV  is  the  volume  of  the  differential  control 
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volume,  if  the  heat  generation  rate,  geometry,  and  heat  convection  rate  were  known, 
the  temperature  gradient  required  could  be  calculated. 


Figure  2.1  Differential  Control  Volume  for  Heat  Flow  Calculations 

2.3.2  Convection.  In  general,  heat  transfer  is  a  case  of  heat  flowing  from  a 
high  temperature  material  to  a  lower  temperature  material.  In  fluid  mechanics  the 
heat  transfer  rate  per  unit  area  (heat  flux)  g”  is  proportional  to  the  difiference  in 
temperature  between  the  fluid  and  the  solid  wall,  as  shown  in  Equation  2.13: 

r  Too)  (2.13) 

The  constant  of  proportionality,  hf,  is  the  heat  transfer  coefficient.  The  challenge  is 
to  accurately  determine  hf  for  the  local  conditions.  The  functions  for  determining 
hf  are  part  of  the  ANSYS®/FLOTRAN™  program  (11),  but  a  check  on  the  output 
for  realism  was  desirable.  For  fully  developed  flow  between  two  flat  plates,  Kays 
and  Crawford  (9:332)  list  the  Nusselt  number  (based  on  hydraulic  diameter),  Nud„ 
{NuDh  =  where  k  is  the  conductivity  of  the  fluid)  as  378  for  a  Reynolds  number 
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of  30,000  and  Prandtl  number  of  0.7.  For  hydraulic  diameter  Dh  {Dh  =  perimeter 
2  •  spacing  for  flat  plates),  of  9  mm  and  thermal  conductivity,  k  =  49.53  milliwatts 
per  meter  Kelvin  (^),  this  would  make  hf  about  2,000  in  fully  developed  flow 
between  flat  plates  4.5  mm  apart.  This  is  somewhat  sensitive  to  temperature,  as  fc  is 
temperature  dependent  for  air,  but  at  temperatures  likely  to  exist  in  a  gas  turbine, 
conductivity,  k,  is  between  25  and  100  so  this  estimate  of  hf  should  not  change 
overly  much. 

A  further  correlation  for  Nusselt  number  is  presented  in  Hill  and  Peterson’s 
book.  Mechanics  and  Thermodynamics  of  Propulsion{20:129).  Again,  this  is  based 
on  fully  developed,  turbulent  flow,  but  in  smooth  tubes.  In  this  correlation  (Equation 
2.14),  the  Nusselt  number  is  proportional  to  the  product  of  the  |  power  of  Reynolds 
number  and  the  cube  root  of  the  Prandtl  number: 

Nud^  =  0.023(Rex),f^(Prf-^^  (2.14) 

Thus,  for  a  given  Reynolds  number  and  Prandtl  number,  the  Nusselt  number  is 
a  constant.  This  does  not,  however,  translate  to  constant  heat  transfer  coefficients. 
As  the  temperature  changes  there  will  be  small  changes  to  the  Reynolds  number  and 
Prandtl  number,  as  well  as  the  thermal  conductivity  of  the  air.  The  equation  can  be 
rearranged  to  give  an  explicit  statement  of  h/  as: 

(Trri  \  -0.2  /  \  -0.67 

(^) 

This  gives  an  additional  estimate  of  the  heat  transfer  coefficients  expected. 
Using  data  for  stagnation  conditions  entering  the  heat  exchanger,  and  a  typical 
velocity:  U  =  29  m/s,  p  is  3.1  •  10“®^,  p  is  6.45^,  Cp  is  lllO^g^,  k  of  0.063 
and  Dh  is  0.009  m.  This  yields  a  hj  of  688  (20:129).  This  is  significantly  lower 

than  the  estimate  from  Kays  and  Crawford,  but  since  these  correlations  are  based 
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on  slightly  different  flows  (smooth  tubes  and  flat  plates),  some  diff'erence  would  be 
expected.  These  two  estimates  did  provide  an  idea  of  the  order  of  magnitude  the 
heat  transfer  coefficients  should  fall  into. 

In  either  case,  most  of  the  flow  for  the  geometries  used  in  this  research  will 
be  thermal  entry  length  flow.  This  should  mean  the  hf  will  be  significantly  higher 
(10-20%  or  more)  at  the  leading  edge  of  the  heat  exchanger  fins,  and  asymptotically 
approach  the  fully  developed  flow  condition  as  the  fin  length  becomes  very  large. 
From  Kays  and  Crawford  (9:340)  we  can  estimate  the  ^  for  the  Nusselt  number  to 
reach  the  fully  developed  value  at  around  35.  Since  the  heat  exchangers  used  in  this 
research  have  a  fin  length  of  0.254  m  and  a  Dh  around  0.009  m,  it  is  easy  to  show 
that  the  ^  will  not  exceed  about  28.  Therefore,  the  flow  will  probably  not  be  fully 
developed,  and  all  of  our  Nud^,  and  thus,  all  hfS  should  exceed  the  fully  developed 
flow  conditions.  There  are  a  number  of  variables  which  affect  the  Nusselt  number 
though.  Increased  surface  roughness  will  tend  to  increase  the  Nusselt  number,  as 
will  the  square  corners  used  on  all  heat  exchanger  fins  in  this  research.  The  heat  flux 
coefficient  data  from  ANSYS®  was  not  expected  to  match  either  of  these  correlations 
exactly,  but  it  was  expected  to  be  of  the  same  order  of  magnitude. 

2.3.3  Pressure  Losses  Due  to  Heating.  There  is  an  inevitable  loss  of  total 
pressure  due  to  heating.  This  can  be  illustrated  through  the  Rayleigh  line  flow,  which 
describes  the  frictionless  flow  of  compressible  fluid  in  a  constant  area,  heated  duct. 
This  can  be  derived  from  a  one-dimensional  momentum  equation  which  reduces  to: 

p  +  pV^  =  constant  (2-16) 

where  p  is  pressure  and  p  is  density.  Since  the  mass  flux  (mass  flow/unit  area),  G, 
is  constant,  pV^  can  be  written  as  GVp,  and  Equation  2.16  can  be  rewritten  as: 

p  -1-  —  =  constant  (2-17) 

P 


2-8 


Therefore,  p  reduces  linearly  with  1/p,  and  since  p  decreases  with  increased 
temperature,  p  must  also  decrease  with  increased  temperature.  To  keep  mass  flow 
constant,  the  speed  must  increase,  which  tends  to  offset  some  of  the  total,  or  stag¬ 
nation,  pressure  loss,  but  this  still  results  in  an  overall  loss  of  stagnation  pressure 
due  to  temperature  increase  (21:205-210;22:249). 

More  complex  formulations  will  give  the  total  pressure  loss  based  on  temper¬ 
ature  change,  inlet  velocity,  etc.  These  are  shown  below. 

Ml  ~  Ml 

Pi  (1+ 7-^1 )  =  P2(1  +  7-^2)  (2.18) 

By  solving  these  equations  for  given  inlet  conditions  and  a  given  outlet  tem¬ 
perature,  the  outlet  Mach  number  can  be  found.  With  the  outlet  Mach  number,  the 
outlet  static  pressure  can  be  caluclated,  and  the  total  pressure  can  be  determined 
using  isentropic  relations.  For  the  flows  considered  in  this  research,  the  temperature 
is  rising  from  622  K  to  approximately  1100  K,  at  inlet  velocities  around  40  m/s  or 
a  Mach  number  of  0.08.  This  combines  to  give  a  total  pressure  loss  of  about  one 
percent,  just  due  to  temperature  increase  (21:205-210). 

2.4  Gas  Turbine  Performance 

2.4.1  On-Design  Gas  Turbine  Performance.  This  section  will  concentrate 
mainly  on  the  effects  of  turbine  inlet  temperature,  Tta,  and  burner  pressure  ratio, 
■Kb,  on  the  thrust  of  a  turbojet  engine.  Using  ONX,  a  software  tool  included  with 
Mattingly,  Heiser  and  Daley’s  Aircraft  Engine  Desipn (4;23),  7x4  and  tt^  were  varied 
to  create  Figure  2.2.  As  shown  in  Figure  2.2,  increases  in  -Kb  and  Tta  both  lead  to 
increased  thrust,  but  there  is  obviously  less  limitation,  and  therefore  more  possible 
gain,  to  increases  in  Tta- 
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— Pressure  Ratio  =  100%  •  J-57  — Pressure  Ratio  =  90% 

. Pressure  Ratio  =  80%  —  Pressure  Ratio  =  70%  --Pressure  Ratio  =  60% 

Pressure  Ratio  =  50%  _ _ _ _ _ 


Figure  2.2  Burner  Pressure  Ratio  and  Turbine  Inlet  Temperature  Effects  on  J-57 
Thrust,  from  ONX 

ONX  estimates  thrust  by  calculating  pressure  and  temperature  at  each  station 
through  the  engine,  using  the  specified  compressor  pressure  ratios,  burner  pressure 
ratios,  turbine  inlet  temperature,  and  the  efficiencies  for  diffusers,  nozzles,  compres¬ 
sors  and  turbines.  The  ONX  program  carries  out  all  of  the  calculations  based  on 
inputs,  calculates  a  work  balance  on  the  high  and  low  pressure  spools  to  determine 
the  pressure  and  temperature  ratios  across  the  high  and  low  pressure  turbines,  deter¬ 
mines  the  fully  expanded  fiow  velocity  at  the  outlet  and  an  exit  area  per  unit  mass 
flow  (23).  These  can  be  combined  with  a  specified  mass  flow  to  calculate  impulse  at 
each  station  through  Equation  2.19: 

It  —  Pi-A-i  -f-  PiV^Ai  (2.19) 
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where  li  is  the  impulse,  P  is  pressure,  A  is  the  area,  p  is  density,  and  V  is  velocity, 
all  at  station  i.  We  can  rearrange  these  terms  by  recognizing  that  m  is  pAV,  to  get: 

U  =  PiAi+miVi  (2.20) 

The  uninstalled  engine  thrust  can  be  calculated  by: 


F  =  Iq  —  lo  —  PoiAg  —  Ao) 


(2.21) 


where  F  is  the  uninstalled  thrust,  subscripts  denote  the  station  number  in  the  engine: 
station  0  is  at  the  beginning  of  a  streamtube  leading  to  the  engine  inlet  (for  subsonic 
conditions)  and  station  9  is  at  the  nozzle  outlet,  for  a  fully  expanded  nozzle.  These 
terms  can  be  further  rearranged,  assuming  constant  mass  flow,  or  very  small  fuel 
flow  relative  to  the  air  flow,  to  get: 


F  =  m 


m 


(2.22) 


Therefore,  since  the  mass  flow  through  the  engine  is  known,  as  are  the  pressures 
and  temperatures,  Mach  number,  velocity  and  area  ratios  can  be  calculated  from 
isentropic  compressible  flow  relationships  (22:110-117).  Since  all  of  the  quantities  on 
the  right  hand  side  of  Equation  2.22  are  known  or  can  be  directly  calculated,  the 
uninstalled  thrust  can  be  easily  calculated.  This  entire  process  was  automated  in 
ONX.  In  this  case,  a  fully  expanded  nozzle  was  assumed,  and,  since  this  is  uninstalled 
thrust,  no  inlet  or  boattail  losses.  This  was  deemed  reasonable  for  the  purposes  of  this 
research,  since  no  airframe  data  was  specified,  and  engine  output  is  being  compared 
to  similar  engines  (4:225;23). 


2.4.2  Off-Design  Gas  Turbine  Performance.  Oif-design  analysis  of  the  op¬ 
eration  of  the  J-57  was  necessary  to  evaluate  the  effect  of  flight  a  higher  altitudes  on 
the  heat  exchanger.  The  program  OFFX,  also  included  with  Mattingly,  Heiser,  and 
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Daley’s  book  (4;23),  was  used  to  accomplish  this.  OFFX  uses  the  engine  specifica¬ 
tions  from  ONX,  including  efficiencies  for  compressors,  turbines,  mechanical  linkages, 
etc;  along  with  the  calculated  high  pressure  turbine  pressure  ratio,  which  is  constant 
when  the  low  pressure  turbine  is  choked,  to  calculate  engine  requirements  away  from 
the  design  point.  By  assuming  that  the  low-pressure  turbine  is  choked  (which  is  true 
over  a  wide  range  of  operating  conditions),  the  high  pressure  turbine  pressure  ratio 
and  corrected  mass  flow  are  held  constant.  OFFX  requires  inputs  for  the  maximum 
allowable  compressor  pressure  ratios,  turbine  inlet  temperature,  total  temperature 
and  pressure  entering  the  burner,  and  high  and  low  pressure  spool  speeds.  OFFX 
can  then  perform  a  work  balance  on  the  high  pressure  spool  to  get  the  high  pressure 
compressor  pressure  ratio,  and  maximize  either  compressor  pressure  ratio  or  turbine 
inlet  temperature,  without  violating  the  other  limitations,  to  calculate  a  thrust,  mass 
flow,  and  other  conditions  through  the  engine,  using  methods  similar  to  ONX.  Using 
Data  from  OFFX,  the  burner  entrance  conditions  at  any  flight  altitude  and  Mach 
number  could  be  found. 

2.5  Triggered  Hafnium  Decay 

The  physics  community  has  an  interest  in  the  4  and  5  quasiparticle  isomers  of 
lutetium  (Lu),  hafnium  (Hf),  and  tantalum  (Ta),  because  they  have  relatively  long 
half-lives,  and  high  excitation  energies.  These  long  half-lives  are  possible  because 
rapid  decay  is  inhibited  by  the  structure  of  the  atomic  nuclei.  Since  rapid  decay  to  a 
lower  energy  state  is  inhibited,  these  isomers  have  a  relatively  long  half-life.  This  is 
interesting  because  this  represents  a  large  amount  of  stored  energy.  In  the  isomer  of 
interest  to  this  research,  ^|Hf,  this  stored  energy  is  approximately  1.3  GJ/g,  around 
one  percent  of  the  energy  of  a  fission  reaction  (2;3). 

Research  shows  that  this  energy  can  be  released  by  bombarding  the  nucleus 
with  X-rays  at  10  to  90  keV.  This  has  the  effect  of  further  raising  the  energy  state  of 
the  nucleus,  and  freeing  the  nucleus  from  the  structural  prohibitions  against  rapid 
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decay.  In  this  manner,  the  stored  energy  can  be  released  a  volumetric  rates  of  up 
to  50  GJ/m^s  (2)(3)(24).  Since  this  reaction  releases  only  photons,  no  particles,  no 
change  to  the  material  properties  of  the  hafnium  was  expected. 


2.6  Radiation  Shielding 

The  first  step  in  determining  the  radiation  shielding  requirements  would  be 
to  determine  the  radiation  output  from  the  source.  In  this  case,  to  run  the  J-57 
engine  full  throttle  at  sea  level,  the  heat  output  to  the  air  must  be  around  48  MW. 
This  should  represent  about  95%  of  the  total  energy  output  of  the  hafnium,  the 
remaining  5%  being  radiation  (24).  Therefore,  the  engine  would  have  around  2.5 
MW  of  radiation  output.  The  equation  for  radiation  dose  is  (25:369): 


b  = 


CE  •  0.0027 

^2 . 4^^ 


(2.23) 


where  CE  is  the  radiation  generation  rate  in  W  and  r  is  the  radius  from  the  source 
in  m.  If  CE  is  2.5  MW  and  r  is  about  3m,  the  dosage,  D  is  about  55  Grays/sec,  or 
approximately  5,500  rem/s.  A  safe  dosage  for  a  human  pilot,  or  similarly  sensitive 
component,  might  be  on  the  order  of  1  ixrem/s  (25:369). 

This  leads  to  the  shielding  requirements  for  the  pilot  or  other  components.  If 
it  is  assumed  that  all  of  the  energy  leaves  the  engine  in  the  form  of  600  keV  gamma 
rays,  deemed  the  most  penetrating  likely  to  escape  the  reaction  process,  and  that 
lead  shielding  is  used  around  sensitive  areas,  by  consulting  tables,  one  can  find  the 
mass  attenuation  (^)(Note  that  within  this  section  only,  n  is  not  viscosity,  but  an 
attenuation  coefficient  for  radiation)  to  be  about  0.13.  Since  the  density  of  lead  is 
about  11.4  this  makes  /x  about  1.5  cm~^.  The  thickness  of  shielding  required 

to  cut  radiation  by  a  given  amount  is  calculated  from  (25:187): 

x=—ln(—]  (2.24) 

/X  \noJ 
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with  fi  in  cm“\  n  the  allowable  dose  behind  the  shielding,  and  rig  the  potential  dose 
without  shielding,  this  gives  the  thickness,  x,  in  cm.  For  this  case,  n  is  10~®  rem/s, 
rio  is  5,500  rem/s,  and  /x  is  1.5  cm"^  making  x  about  15  cm,  or  a  shield  load  of  1,700 
kg/rri^,  or  350  Ib/ff  of  shield  required.  This  should  not  be  taken  as  an  absolute 
or  very  accurate  number.  This  is  a  very  low  order  estimate,  which  mainly  goes  to 
show  the  need  for  a  significant  mass  of  radiation  shielding.  The  estimate  above  could 
be  off  by  a  significant  factor,  depending  on  the  exact  spectrum  of  radiation  to  be 
attenuated,  the  geometry  involved,  the  exact  choice  of  materials,  and  the  overall 
design  of  the  shield.  Therefore,  all  that  can  be  said  with  accuracy  at  this  point  is 
that  the  weight  of  the  radiation  shielding  will  be  significant  (hundreds  or  thousands 
of  pounds),  and  that  it  will  be  highly  airframe  dependent  (25:187). 
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III.  Computational  Setup 

This  chapter  contains  a  description  of  the  computational  hardware,  software, 
and  methods  used  in  this  research. 

3. 1  Hardware 

All  computations  for  this  research  were  accomplished  on  a  Dell®  Dimension 
XPS  T700R  personal  computer.  This  computer  was  equipped  with  a  700  MHz 
Pentium  III®  processor,  384  Mb  of  RAM,  and  a  19”  monitor. 

3.2  Software 

This  research  used  the  commercial  software  package,  ANSYS®  5.6.1  (11).  This 
package  includes  the  FLOTRAN®  computational  fluid  dynamics  program.  ANSYS® 
comes  equipped  with  its  own  grid  generation  and  is  capable  of  analyzing  many  types 
of  flows,  from  simple  adiabatic,  incompressible,  laminar  flows,  to  turbulent,  com¬ 
pressible  flows  with  heat  transfer.  Multiple  solver  options,  advection  options,  tur¬ 
bulence  models  and  wall  treatments  are  available.  Boundary  conditions  on  velocity, 
pressure,  temperature,  heat  flux,  and  heat  generation  are  available  and  were  used  in 
this  research. 

3.3  Software  Inputs 

This  section  contains  information  about  the  various  software  settings  used  in 
this  research.  For  an  example  of  the  steps  and  inputs  required  to  create  a  heat 
exchanger  model  like  those  used  in  this  research,  see  Appendix  B. 

3.3.1  Solver  Settings.  For  all  runs  the  following  settings  were  used.  All 
runs  were  transient  solutions  (run  to  steady  state),  compressible,  thermal,  and  tur¬ 
bulent.  The  turbulence  model  used  in  all  cases  was  the  k  —  e  from  Shi,  Zhu  and 
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Lumley  (19),  with  the  default  constants.  The  solver  was  set  for  5,000  time-steps, 
with  20  global  iterations  per  time  step,  and  the  ANSYS®  (11)  default  time-step  ter¬ 
mination  criteria.  All  fluid  properties  were  set  to  calculate  in  SI  units,  and  reference 
conditions  were  set  to  match  the  J-57  compressor  outlet  conditions  (shown  later  in 
Table  4.1).  Stability,  relaxation,  and  results  capping  values  were  left  at  the  default 
values.  The  various  equations  used  the  following  solvers:  X- velocity,  Y-velocity,  and 
temperature  used  the  tri-diagonal  matrix  algorithm  solver  (TDMA);  pressure  used 
the  preconditioned  conjugate  residual  solver  (PCRS);  and  turbulent  kinetic  energy 
and  dissipation  equations  used  the  preconditioned  generalized  minimum  residual 
method  (PGMR)  (26:7.39,40).  The  advection  options  selected  were:  for  momen¬ 
tum,  compressible  pressure  and  energy  equations,  the  monotone  streamline  upwind 
method;  and  for  turbulent  equations,  the  streamline  upwind  Petrov/Galerkin  method 
(26:7.21-24).  The  various  settings  were  selected  for  the  stability  they  demonstrated 
in  the  software  verification  process. 

3.3.2  Geometry.  Three  basic  geometries  of  heat  exchanger  were  explored 
in  this  research:  Concentric  annular  tubes,  radial  trapezoidal  fins,  and  a  dual  an¬ 
nulus  of  rectangular  fins.  These  were  selected  as  simple  geometries  which  should 
be  manufactureable.  Simplicity  of  form  may  be  critical,  since  all  of  the  heated  sur¬ 
faces  would  be  made  of  the  hafnium  isomer  and  its  alloyed  materials,  and  thus, 
the  heat  exchanger  would  generate  heat  within  the  solid  material  of  the  heat  ex¬ 
changer.  Given  the  extraordinary  strength  of  hafnium  (yield  strength  about  50% 
higher  than  titanium,  ulitmate  strength  over  100%  higher  (27:665))  machining  may 
be  quite  difficult,  so  it  would  probably  be  helpful  to  keep  manufacturing  steps  to  a 
minimum.  In  addition,  the  flow  through  all  of  these  geometries  can  be  reasonably 
approximated  by  a  series  of  two-dimensional  flow  fields.  This  was  important,  as  a 
single  three-dimensional  run  could  have  taken  several  weeks  on  a  PC,  or  consumed 
a  considerable  amount  of  (very  expensive)  time  on  a  faster  computer  system. 
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3.3.3  Concentric  Tubes.  With  the  concentric  tubes,  the  flow  is  through  an 
annulus  made  up  of  a  number  of  concentric  tubes,  surrounding  the  drive  spools  of  the 
engine.  A  schematic  view  of  the  cross  section  at  section  A-A  of  Figure  1.2  is  shown  in 
Figure  3.1.  The  flow  between  any  two  of  the  concentric  tubes  can  be  approximated 
as  the  flow  between  two  flat  plates,  provided  the  ratio  of  the  radii  of  the  cylinders, 
is  close  to  unity.  Since  the  spacing  used  is  small,  and  there  is  a  relatively  large 
minimum  inner  radius  for  the  heat  exchanger,  this  value,  at  a  minimum  is  about  0.96 
(9:327).  Therefore,  one  can  consider  this  a  two  dimensional  flow  with  a  streamwise, 
or  axial  direction  and  a  radial  direction.  For  this  flow,  two  heat  exchanger  envelopes 
were  used.  The  first,  used  for  many  runs  to  reduce  computation  time,  used  a  25.4 
mm  high  inlet,  expanding  in  76.2  mm  to  a  50.8  mm  high  straight  duct,  255  mm 
long,  narrowing  back  to  a  25.4  mm  outlet,  and  attached  to  a  long  outlet  duct  (see 
Figures  3.2,  3.3  and  3.4).  The  small  size  of  the  geometry  allowed  a  smaller  number  of 
elements,  and  a  vastly  reduced  computation  time  to  evaluate  trends.  This  geometry 
did  not  correlate  to  anything  which  could  be  installed  in  a  J-57  engine,  and  this  was 
intentional.  It  was  used  solely  as  a  tool  for  evaluating  some  trends  and  ideas  in  a 
timely  manner.  In  this  geometry,  shown  in  Figures  3.2,  3.3  and  3.4  as  the  view  of 
Section  A-A  of  Figure  3.1,  tubes  from  0.5  mm  thick  to  2.0  mm  thick  were  checked 
with  about  4.5  mm  spacing  between  them.  With  0.5  mm  tubes,  both  full  length 
tubes,  254  mm  long,  and  staggered  tubes,  10  rows  of  25.4  mm  tubes  offset  2.5  mm 
from  each  other  were  used. 

The  second  geometry  was  generated  to  more  closely  match  the  actual  size  of 
the  J-57  combustion  chamber,  and  Figure  3.5  shows  the  view  from  Section  A-A  of 
Figure  3.1  for  this  geometry.  This  layout  uses  a  50.8  mm  high  inlet,  expanding  in 
254  mm  to  a  206  mm  high  straight  duct,  254  mm  long,  which  contracts  in  254  mm 
to  a  76.2  mm  high  outlet,  connected  to  a  long  outlet  duct.  This  used  2.0  mm  thick 
tubes,  spaced  4.5  mm  apart  for  the  length  to  the  206  mm  high  duct  section,  to  verify 
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Flow  Direction  ^ 


Figure  3.3  Model  Geometry,  2.0  mm  Thick  Tubes,  Compact  Concentric  Geometry 
(Cross-Section  A- A  of  Figure  3.1) 


Figure  3.4  Model  Geometry,  0.5  mm  Thick  Staggered  Tubes,  Compact  Concentric 
Geometry  (Cross-Section  A-A  of  Figure  3.1) 


the  applicability  of  results  from  the  smaller  geometry,  and  evaluate  altitude  effects 
on  the  heat  exchanger. 


Figure  3.5  Model  Geometry,  2.0  mm  Thick  Tubes,  Full-Scale  Concentric  Geometry 
(Cross-Section  A-A  of  Figure  3.1) 

3.5.4  Trapezoidal  Fins.  The  use  of  straight  fins,  for  easier  manufacturing, 
presented  two  options:  either  the  thickness  of  the  fins,  or  the  space  between  the 
fins,  must  change  with  radius.  For  the  trapezoidal  fins,  the  thickness  of  the  fins 
is  changed  with  the  spacing  held  constant.  This  is  shown  schematically  in  Figure 
3.6  as  the  view  from  Section  A-A  of  Figure  1.2,  and  an  expanded  view  in  Figure 
3.7.  This  was  treated  as  several  cases  of  flow  between  flat  plates,  at  the  locations 
indicated  in  Figure  3.7(sections  A-A,  B-B,  and  C-C).  The  geometry  of  a  flat  plate 
case  is  shown  in  Figure  3.8  and  all  radial  fin  geometries  are  similar.  The  dimensions 
of  the  geometry  for  each  section  are  found  in  Table  3.1  (see  Figure  3.8  for  locations). 

3.3.5  Rectangular  Fins.  With  the  rectangular  fins,  two  annuli  of  straight, 
constant  thickness  fins  were  used.  The  view  of  Section  A-A  of  Figure  1.2  is  shown 
in  Figure  3.9,  and  a  close  view  of  a  pair  of  rectangular  fins  is  shown  in  Figure  3.10. 
Two  annuli  were  used  to  minimize  the  variation  in  spacing  between  the  fins.  For 
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Table  3.1  Dimensions  of  Geometries  Used  in  Trapezoidal  Fin  Computations  (see 

_ Figure  3.8  for  Locations  of  Dimensions) _ _ _ 

I  T,nefl.tinn  I  Inlet  H  I  Inlet  L  I  Fin  Soacing  I  Fin  L  Outlet  H  I  Outlet  H 


Inner  Rad. 

5.0  mm 

25.4  mm 

4.5  mm 

254  mm 

5.0  mm 

220.6  mm 

Mid.  Rad. 

6.5  mm 

25.4  mm 

4.5  mm 

254  mm 

6.5  mm 

220.6  mm 

Outer  Rad. 

8.3  mm 

25.4  mm 

4.5  mm 

254  mm 

8.3  mm 

220.6  mm 

Isomer  Fins 


Section  A-A:  Inner  Radius  of  Annulus 
Section  B-B:  Mid-Radius  of  Annulus 
Section  C-C:  Outer  Radius  of  Annulus 


Flow  is  into  Page 


Figure  3.7  Schematic  View  of  Two  Trapezoidal  Fins 


Figure  3.8  Example  of  Radial  Fin  Computational  Geometry 
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each  annulus,  runs  were  done  at  the  inner  radius  (section  A-A  of  Figure  3.10),  outer 
radius  (section  B-B  of  Figure  3.10),  and  center  (section  C-C  of  Figure  3.10)  of  each 
annulus,  to  establish  a  variation  across  the  annulus.  Construction  of  each  model  was 
similar  to  the  trapezoidal  models  (see  Figure  3.8),  except  that  ledges  were  all  1.0 
mm  high,  for  a  2.0  mm  thick  fin,  and  spacing  was  varied.  Dimensions  are  found  in 
Table  3.2,  for  each  of  the  rectangular  fin  computational  models. 

Table  3.2  Dimensions  of  Geometries  Used  in  Rectangular  Fin  Computations  (see 


Figure  3.8  for  Locations  of  Dimensions) 


Location 

Inlet  H 

Inlet  L 

Fin  Spacing 

Fin  L 

Outlet  H 

Outlet  H 

Inner  Annulus 

Inner  Rad. 
Mid.  Rad. 
Outer  Rad. 

5.55  mm 
6.5  mm 
7.45  mm 

25.4  mm 
25.4  mm 
25.4  mm 

3.55  mm 
4.5  mm 
5.45  mm 

254  mm 
254  mm 
254  mm 

5.55  mm 
6.5  mm 
7.45  mm 

220.6  mm 
220.6  mm 
220.6  mm 

Outer  Annulus 

Inner  Rad. 
Mid.  Rad. 
Outer  Rad. 

5.77  mm 
6.5  mm 
7.22  mm 

25.4  mm 
25.4  mm 
25.4  mm 

254  mm 
254  mm 
254  mm 

5.77  mm 
6.5  mm 
7.22  mm 

220.6  mm 
220.6  mm 
220.6  mm 

Figure  3.9  Schematic  View  of  Rectangular  Radial  Fin  Geometry 
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Isomer  Fins 


Section  B-B:  Mid-Radius  of  Annulus 
Section  C-C:  Outer  Radius  of  Annulus 

Flow  is  into  Page 


Figure  3.10  Schematic  View  of  Two  Rectangular  Fins 
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3.3. 6  Radial  Fin  Inlet  and  Outlet.  These  two  sections  were  done  to  provide 
estimates  of  the  pressure  losses  heading  into  and  out  of  the  radial  fin  geometries.  The 
inlet  and  outlet  were  created  to  match  the  J-57  geometry.  The  inlet  was  254  mm 
long,  expanding  the  fiow  from  a  50.89  mm  inlet  height  to  a  206  mm  outlet  height 
(see  Figure  3.11).  The  outlet  section  was  254  mm  long,  compressing  the  flow  from 
the  206  mm  height  to  a  76.2  mm  outlet  (see  Figure  3.12).  Each  of  these  had  a 
746  mm  straight  duct  attached  to  the  outlet  to  allow  the  outlet  pressure  boundary 
condition  to  be  separated  from  the  area  of  interest.  Flow  data  was  taken  at  the  inlet 
and  across  the  exit  of  the  254  mm  long  area  of  interest. 


Figure  3.11  Radial  Geometry  Inlet 

3.3.7  Grid  Generation.  All  computational  domains  were  created  using 
the  built-in  tools  of  ANSYS®/FLOTRAN^^  (11).  The  meshes  were  constructed 
using  unstructured  triangular  elements.  These  were  selected  due  to  the  irregular 
geometries  of  many  domains,  and  a  restriction  in  the  software  against  meshing  some 
irregular  shapes  with  quadrilaterals.  There  are  potential  drawbacks  to  this  type  of 
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Figure  3.12  Radial  Geometry  Outlet 


mesh,  but  to  minimize  these,  element  sizes  were  kept  small  in  the  areas  of  interest 
to  keep  values  in  the  10  -  5,000  range  recommended  by  ANSYS® 

(28).  Outside  the  areas  of  interest,  in  the  exit  ducts  from  the  heat  exchanger  sections, 
elements  were  made  as  large  as  possible  to  decrease  computational  time.  Since  these 
ducts  were  only  present  to  move  boundary  conditions  well  away  from  the  points  were 
data  was  taken,  this  was  judged  to  have  few  drawbacks. 

3.3.8  Reference  Conditions.  All  reference  conditions  were  calculated  from 
the  compressor  output  conditions  for  the  J-57  engine,  using  OFFX  for  operating 
conditions  other  than  static,  sea  level.  Reference  pressure  was  1.15  MPa  for  static  sea 
level  runs  (167  psia),  reference  temperature  was  622  K  (1120  "R)  (29:474).  All  fluid 
properties  (viscosity,  density,  conductivity)  were  calculated  based  on  the  computed 
static  pressure  and  temperature  at  each  location.  However,  ANSYS®  (11)  did  not 
allow  variable  speciflc  heats  with  a  compressible  flow.  Therefore,  a  specific  heat  of 
1004  J/kgK  was  used,  and  the  data  was  corrected  at  the  end  to  account  for  the 
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specific  heat  variation.  The  correction  process  can  be  found  in  Appendix  A.  Using 
OFFX,  combustion  chamber  entrance  conditions  were  calculated  at  Mach  0.8  for 
altitudes  of  15,000  ft,  30,000  ft  and  45,000  ft. 

3.3.9  Boundary  Conditions.  For  all  runs,  similar  boundary  conditions  were 
applied.  First,  all  solid  boundaries  were  represented  by  conditions  of  zero  velocity 
in  both  the  x  and  y  directions.  Inlet  velocity  was  entirely  in  the  x  direction.  The 
compact  concentric  geometries  were  run  with  inlet  velocities  of  75  m/s,  and  in  the 
full  scale  concentric  fin  case,  the  inlet  velocity  was  determined  from  the  J-57  mass 
flow  and  the  inlet  annulus  area  (112  m/s  at  sea  level,  104  m/s  at  15,000  ft,  99 
m/s  at  30,000  ft,  and  95  m/s  at  45,000  ft).  For  the  radial  fin  geometries,  the  inlet 
velocities  were  run  at  two  levels:  the  first,  at  41  m/s  was  to  closely  match  the 
compact  concentric  geometry  results  for  comparison;  the  second,  at  29  m/s  was  to 
match  the  one-dimensional  expansion  from  the  J-57  burner  inlet  to  the  nominal  heat 
exchange  height  (a  1:4  ratio).  The  second  velocity  should  closely  match  the  actual 
velocities  in  the  J-57  burner  which  should  25-30  m/s  at  sea  level  static  operation.  At 
the  outlet  from  the  models,  pressure  was  set  to  zero,  relative  to  the  pressure  offset 
input. 


3.3.10  Thermal  Boundary  Conditions.  To  establish  the  proper  heat  trans¬ 
fer  boundary  condition  for  the  flow,  three  possible  cases  were  evaluated:  constant 
wall  temperature,  constant  heat  flux  (at  the  fin  or  tube  surface),  and  constant  volu¬ 
metric  heat  generation  (within  the  fin  or  tube  material).  Since  the  hafnium  isomer 
was  assumed  to  make  up  the  entirety  of  the  heat  exchanger  fins  or  tubes,  the  heat 
would  be  generated  within  the  body  of  the  fins  or  tube  themselves.  The  fins  or  tubes 
were  assumed  to  be  manufactured  from  hafnium  because  hafnium  is  a  structurally 
strong  metal  (the  yield  strength  of  hafnium  is  50%  greater  than  that  of  titanium) 
(27:665),  and  the  surface  temperatures  involved  would  melt  most  other  materials. 
The  constant  heat  generation  was  assumed  to  be  the  most  accurate  model  of  the 
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isomer  decay  process,  and  the  best  heat  transfer  boundary  condition.  However,  it  is 
also  the  most  computationally  intensive,  as  it  requires  the  temperature  to  stabilize 
in  the  material,  which  takes  a  significant  amount  of  time  (An  early  estimate  with 
small  models  indicated  that  the  run  time  for  a  constant  heat  generation  case  would 
be  approximately  one  month).  Therefore,  it  was  decide  to  select  one  of  the  other 
boundary  conditions,  which  do  not  require  convergence  to  their  own  steady  state 
before  the  fiow  can  begin  to  converge.  To  decide  between  constant  wall  temperature 
or  constant  heat  flux,  it  was  decided  to  go  with  the  condition  whose  temperature 
gradient  most  nearly  matched  the  constant  heat  generation  condition.  The  variation 
of  temperature  distributions  is  shown  in  Figure  3.13. 


—  Constant  Heat  Generation  —  Constant  Heat  Flux 
. Constant  Wall  Temperature _ 


Figure  3.13  Variation  of  Wall  Temperature  along  Fin/Tube  Surface 

The  constant  wall  temperature,  of  course  has  no  variation  in  temperature  along 
the  plate.  The  constant  heat  flux  shows  a  variation  of  about  660  K/m,  and  the 
constant  heat  generation  shows  a  variation  of  about  120  K/m.  The  mechanism  for 
the  difference  between  constant  heat  flux  and  constant  wall  temperature  was  the 
conduction  through  the  material  of  the  fin.  The  trailing  edge  of  the  fin  generates 
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heat  which  cannot  be  removed  by  convection,  since  the  heat  transfer  coefficients  are 
lower,  so  this  heat  was  conducted  toward  the  leading  edge,  where  the  heat  transfer 
coefficients  are  higher,  and  the  heat  can  be  convected  away.  For  all  computations 
in  this  research,  a  constant  wall  temperature  of  2400  K  was  applied  to  the  surfaces 
of  the  hafnium  isomer  fins  or  tubes,  and  all  other  solid  surfaces  were  considered 
insulated  (zero  heat  flux  boundary  condition  applied).  This  temperature  (2400  K) 
was  chosen  as  a  reasonable  margin  below  the  melting  point  of  hafnium  (around  2500 
K)  (27:665). 

It  is  worth  noting  that,  where  there  are  significant  changes  in  fin  thickness, 
this  may  be  a  very  poor  assumption.  If  an  energy  balance  is  performed  on  a  small 
control  volume  of  fin,  as  shown  in  Figure  2.1,  it  can  be  shown  that,  if  the  net  increase 
in  stored  energy  within  the  control  volume  is  zero  (steady  state  operation): 


^generated  d"  Qcond.in  Qcond.out  Qconv.  0 

where  these  quantities  are  given  by: 

Qgenerated 
Qcond,in 

^COTld.Q'xif, 

Qconv. 

and  HGR  is  the  volumetric  heat  generation  rate.  If  HGR  is  taken  to  be  1  ■  10®  ^ 
(the  approximate  level  required  for  J-57  operation  at  sea  level  with  the  trapezoidal 
fin  geometry).  A;  =  22  and  h/  =  700  (a  reasonable  average  value  for  the 
observed  data)  it  can  be  quickly  shown  that  very  large  magnitudes  of  ^  are  required 
at  both  extremes  of  the  trapezoidal  fins  to  balance  the  energy  in  the  control  volume. 
At  the  thick  end  (t  =  3.8  mm),  ^  must  be  large  to  conduct  heat  away  from  the 


{HGR)  ■  '^^dr 


t 


k 


dr 


(t  T  dt) 
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control  volume,  while  at  the  thin  end  (t  =  0.5  mm),  §  must  be  very  large  to 
conduct  heat  into  the  control  volume.  In  either  case,  to  assume  constant  volumetric 
heat  generation  and  constant  wall  temperatures  in  the  radial  direction,  the  fins  must 
be  nearly  constant  thickness,  or  they  must  be  hollow,  or  filled  with  neutral  material, 
such  that  the  volume  distribution  of  isomer  is  nearly  constant  at  any  radius  along 
the  fin.  If  solid  fins,  varying  in  thickness  from  0.5  mm  to  3.8  mm  were  used,  the 
temperatures  could  vary  hundreds  to  thousands  of  degrees  from  the  inner  radius  to 
the  outer  radius.  In  the  constant  thickness  fins,  there  must  be  some  temperature 
gradient  (^)  to  account  for  heat  transfer  from  the  inner  radius  (where  hf  would 
be  low)  to  the  outer  radius  (where  hf  would  be  higher),  but  these  would  be  much 
smaller  (on  the  order  of  100  K  across  the  radial  span  of  each  fin  (103  mm),  and 
therefore  a  maximum  temperature  difference  on  the  order  of  50  K  at  any  point),  and 
negligible  in  light  of  the  normal  accuracy  of  heat  transfer  calculations. 

Therefore,  while  constant  temperature  in  both  the  axial  and  radial  directions 
is  a  reasonable  assumption  for  constant  (or  nearly  constant)  thickness  fins  or  tubes, 
if  the  thickness  changes  significantly  (as  in  the  constant  spacing  radial  fin  geometry), 
this  may  be  very  misleading.  The  exception  to  this  would  be  a  tapered  fin  which  is 
manufactured  hollow  (forming  a  V  in  cross  section),  so  that  the  total  thickness  of 
material  is  nearly  constant  in  the  radial  direction.  In  this  research,  to  assume  solid 
trapezoidal  fins  (of  the  dimensions  specified)  would  probably  require  fin  temperatures 
in  excess  of  the  melting  point  of  hafnium  to  maintain  the  overall  levels  of  heat  transfer 
obtained  here.  The  data  presented  in  this  research  for  that  geometry  would  only  be 
possible  with  fins  made  up  of  two  thin  plates  forming  the  trapezoidal  fin  (possibly 
filled  with  some  material  which  would  not  contribute  to  the  heat  generation,  but 
which  would  contribute  structural  support). 
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3-4  Code  Validation 

In  preparation  for  conducting  this  research,  ANSYS®/FLOTRAN®  (11)  was 
used  to  replicate  the  results  of  research  conducted  by  Jarcy  (30)  in  evaluating  the 
flow  fleld  behind  a  heated  cylinder  in  crossflow.  Using  identical  geometry  and  test 
conditions,  a  very  similar  total  pressure  profile  was  obtained  behind  the  cylinder. 
The  profile  obtained  was  judged  even  better  when  the  effects  of  the  relatively  large 
pitot  probe  used  by  Jarcy  to  obtain  his  data  were  considered.  With  the  inevitable 
smearing  of  sharp  pressure  gradients  induced  by  the  pitot  probe,  the  experimental 
data  and  computational  data  are  very  close.  All  of  the  settings  feeding  into  the 
turbulence  model  were  maintained  throughout  the  research  (30;31). 


Figure  3.14  Comparison  of  Experimental  and  Computational  Dynamic  Pressure 
Distributions  Behind  a  Cylinder  in  Crossflow 
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Table  3.3  Comparison  of  Temperature  Distributions  on  a  Cylinder  in  Crossflow 


Location  in 
Degrees  from 
Stagnation  Point 

Experimental 
Temperature  (®F) 

Computational 
Temperature  (°F) 

0 

553 

561 

90 

734 

651 

180 

739 

728 
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IV.  Results  and  Discussion 

This  chapter  contains  the  results  of  the  computational  models  in  the  present 
work.  The  estimated  performance  of  the  various  heat  exchanger  geometries  m  a  J-57 
turbojet  is  also  provided. 

4.1  Pressure  Ratios  and  Total  Temperature  Increases 

All  data  is  provided  in  terms  of  average  total  pressure  ratio  (expressed  in 
a  percentage,  outlet  over  inlet)  and  an  average  total  temperature.  For  the  reduced 
size  geometries  for  concentric  tubes,  these  are  algebraic  averages,  while  for  all  others, 
an  area-weighted  average  is  used  This  should  provide  the  best 

estimate  of  the  turbine  inlet  conditions  for  the  operating  gas  turbine  engine.  Note 
that  all  temperatures  in  charts  are  uncorrected,  but  corrected  temperatures  have 
been  used  in  all  performance  estimates.  Overall  results  (corrected  temperatures  and 
pressure  ratios)  can  be  found  in  Table  4.6. 

4.1.1  Concentric  Tube  Geometry.  The  first  section  run  used  the  compact 
geometry  for  0.5  mm  thick  concentric  tubes.  Again,  this  geometry  was  not  intended 
to  correlate  to  anything  installed  in  the  J-57  turbojet,  it  was  solely  intended  to 
examine  trends  and  options  in  a  timely  manner  (the  reduced  size  led  to  reduced 
computation  time).  The  outlet  conditions  for  this  are  shown  in  Figure  4.1,  as  they 
vary  from  the  inner  (lower)  wall  of  the  geometry  shown  in  Figure  3.2  (location  C) 
to  the  upper  wall  (location  D).  As  one  would  expect,  the  total  pressure  is  nearly 
constant  across  the  outlet,  while  the  total  temperature  has  not  yet  mixed  enough  to 
reach  a  near  uniform  condition.  The  (algebraic)  average  outlet  conditions  for  this 
geometry  are:  Tt4  =  1090  K  (1962‘’R),  and  tt;,  =  96.76%.  This  Tn  corrects  to  1027 
K  (1848‘’R)  (see  Table  4.6  for  comparison  to  other  geometries). 
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An  item  of  interest  was  to  evaluate  the  heat  transfer  coefficients  for  reason¬ 
ableness.  Figure  4.2  is  a  plot  of  the  computed  heat  flux  and  heat  transfer  coefficient 
variation  along  the  center  (fifth  from  the  bottom)  tube/fin  of  the  geometry  shown 
in  Figure  3.2,  starting  at  the  left  (location  A)  and  going  along  the  upper  surface  of 
the  tube/fin  to  location  B.  As  shown  in  Figure  4.2,  the  heat  transfer  coefficient  has 
a  high  of  about  1,000  Wjm^K  and  a  low  of  about  700  This  is  well  short 

of  the  value  of  about  2,000  W/m^K  for  fully  developed  turbulent  flow  between  two 
flat  plates  predicted  by  Kays  and  Crawford  (9:326).  However,  it  is  very  close  to  the 
value  of  688  W/rn^K  calculated  from  Hill  and  Peterson  (  See  Section  2.3.2  of  this 
document)(20:129).  The  leading  edge  hf  is  about  30%  higher  than  the  trailing  edge, 
and  the  data  shows  overall  excellent  agreement  with  the  expected  results.  Therefore, 
the  ANSYS®  (11)  data  was  taken  as  a  reasonable  approximation  of  reality.  Other 
concentric  tube  geometries  showed  behavior  similar  to  Figure  4.2. 


Figure  4.1  Radial  Variation  of  Total  Temperature  and  Total  Pressure  Across  Exit 
Plane  (Location  C  to  Location  D  in  Figure  3.2),  Compact  Concentric 
Tube  Heat  Exchanger,  0.5  mm  Thick  Tubes 


Figure  4.2  Axial  Variation  of  Heat  Flux  and  Heat  Transfer  Coefficient  on  Center 
Tube  (Location  A  to  Location  B  in  Figure  3.2),  Compact  Concentric 
Tube  Heat  Exchanger,  0.5  mm  Thick  Tubes 

Two  additional  configurations  were  run  with  the  compact  concentric  tube  ge¬ 
ometry  (see  Figures  3.1,  3.2,  and  3.3  for  illustrations).  The  first  used  fins  2.0  mm 
thick  instead  of  0.5  mm  thick  (see  Figure  3.3).  This  was  done  to  evaluate  the  pressure 
and  temperature  effects  of  thicker  tubes,  as  might  thicker  tubes  might  be  required 
for  structural  integrity  or  ease  of  manufacture.  With  the  total  heat  exchanger  height 
the  same,  and  maintaining  constant  spacing  between  the  tubes,  this  resulted  in  fewer 
tubes.  Results  are  shown  in  Figure  4.3,  as  the  total  temperature  and  total  pressure 
distributions  across  the  exit  of  the  heat  exchanger  from  the  lower/inner  wall  (loca¬ 
tion  A  in  Figure  3.3)  to  the  upper  wall  (location  B).  The  average  conditions  were: 
Tta  =  1042  K  (ISTfi^R),  and  itb  =  97.01%.  This  corrects  to  986  K  or  1774°R  (see 
Table  4.6  for  comparison  with  other  geometries).  The  Tta  seems  low,  but  some  of 
this  is  due  to  the  space  between  the  inner  wall  and  the  first  tube  being  somewhat 
larger.  If  the  outer  half  of  the  data  points  were  used,  the  average  Tta  would  be  1092 
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K  (1966°R) (which  corrected  to  1028  K  or  ISSPR).  The  somewhat  lower  pressure 
losses  can  be  attributed  to  the  lower  temperature  rise,  combined  with  lower  viscous 
losses  due  to  less  exposed  surface  area. 


Distance  From  Inner  Wall  (m) 

I - Total  Pressure  —Total  Temperatur^ 

Figure  4.3  Radial  Variation  of  Total  Temperature  and  Total  Pressure  Across  Exit 
Plane  (Location  A  to  Location  B,  in  Figure  3.3),  Compact  Concentric 
Tube  Heat  Exchanger,  2.0  mm  Thick  Tubes 

The  last  configuration  run  with  the  compact  concentric  tubes  uses  staggered 
rows  of  tubes,  each  25.4  mm  long,  offset  by  2.5  mm  (this  geometry  is  shown  in  Figure 
3.4).  This  was  done  to  evaluate  the  possibility  of  significantly  increased  heat  transfer. 
Higher  heat  transfer  was  expected  because  almost  the  entire  surface  area  of  the  heat 
exchanger  would  be  in  the  very  beginning  of  the  thermal  entry  length,  since  each 
row  of  tubes  would  start  a  new  entry  length.  Therefore,  the  heat  exchanger  would 
have  the  highest  possible  heat  transfer  coefficients  throughout.  Results  are  shown 
in  Figure  4.4  (total  temperature  and  total  pressure  distributions  from  location  A 
in  Figure  3.4  to  location  B),  and  the  average  outlet  conditions  are:  Tta  =  1108  K 
(1994°R)  and  ttj  =  95.0%.  The  corrected  Tta  is  1042  K,  or  1875“R.  There  is  a  small 
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improvement  in  heat  transfer,  but  the  relatively  large  increase  in  pressure  loss  and 
increased  complexity  would  likely  offset  the  small  possible  performance  gains  (refer 
to  Table  4.6  for  performance  comparisons). 


Figure  4.4  Radial  Variation  of  Total  Temperature  and  Total  Pressure  Across  Exit 
Plane  (Location  A  to  Location  B  in  Figure  3.4),  Compact  Concentric 
Tube  Heat  Exchanger,  0.5  mm  Thick  Staggered  Tubes 

A  final  geometry  was  evaluated  with  the  concentric  tubes.  This  was  a  full  scale 
model,  which  would  simulate  an  annular  heat  exchanger  fully  replacing  the  burner 
section  of  a  J-57  (see  Figures  3.1  and  3.5.  The  overall  length  is  76.2  cm  (30  inches), 
and  the  heat  exchanger  has  a  maximum  height  of  20.6  cm  (8  inches).  The  tubes  are 
2  mm  thick,  4.5  mm  apart,  and  25.4  cm  long  (10  inches).  Higher  total  temperatures 
were  expected  out  of  this  geometry  for  two  reasons:  First,  the  cooling  passages  at 
the  upper  and  lower  wall  make  up  far  less  of  the  total  flow;  Second,  the  flow  entering 
the  fins  should  be  a  lower  velocity  than  evaluated  in  the  compact  geometry.  Even 
though  the  lower  velocity  would  indicate  a  lower  Reynolds  number,  and  therefore 
lower  heat  transfer  coefficients,  we  would  expect  the  outlet  temperature  to  increase. 
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since  the  mass  flow  through  any  given  channel  decreases  faster  than  the  heat  transfer 
coefficient  (hf  <x  Re%l  while  m  a  Since  the  total  heat  added  to  the  flow 

does  not  decrease  as  fast  as  the  mass  flow,  the  temperatures  must  increase.  Higher 
pressure  losses  were  also  expected,  due  to  larger  temperature  increases,  and  higher 
diffusion/compression  ratios  for  the  inlet  and  outlet  manifolds,  although  this  would 
be  balanced  somewhat  by  the  lower  viscous  losses  due  to  the  lower  velocities  on  the 
tubes.  Overall  distributions  for  outlet  total  temperature  and  pressure  are  shown 
in  Figure  4.5,  from  the  inner  radius  (location  A  in  Figure  3.5)  to  the  outer  radius 
(location  B  in  Figure  3.5).  The  average  (area  weighted)  Tta  was  1236  K  (2224'’R) 
and  TTft  was  94.5%.  This  led  to  a  corrected  Tta  of  1150  K  (2070°R)  (see  Table  4.6  for 
performance  comparisons  to  other  geometries). 


Figure  4.5  Radial  Variation  of  Total  Temperature  and  Total  Pressure  Across  Exit 
Plane  (Location  A  to  Location  B  in  Figure  3.5),  Full-Scale  Concentric 
Tube  Heat  Exchanger,  2.0  mm  Thick  Tubes 

The  full-scale  concentric  tube  geometry  was  also  evaluated  a  three  altitudes 
to  establish  the  ability  of  the  heat  exchanger  to  function  over  a  reasonable  flight 
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envelope.  Data  from  OFFX  (23)  was  used  to  generate  new  inlet  conditions  for  the 
three  different  cases:  15,000  ft  at  Mach  0.8;  30,000  ft  at  Mach  0.8;  and  45,000  ft  at 
Mach  0.8.  These  are  found  in  Table  4.1. 


Table  4.1  J-57  Burner  or  Heat  Exchanger  Inlet  Conditions  at  Typical  Flight  Alti- 


tiides  and  Mach  Numbers 


Altitude 

Mach  Number 

Total  Temperature 

Total  Pressure 

Inlet  Velocity 

Oft 

0 

622  K 

1.15  MPa 

112  m/s 

15,000  ft 

0.8 

625  K 

0.97  MPa 

104  m/s 

30,000  ft 

0.8 

566  K 

0.53  MPa 

99  m/s 

45,000  ft 

0.8 

536  K 

0.26  MPa 

95  m/s 

Each  set  of  inlet  conditions  was  evaluated  for  its  effect  on  (corrected)  Tt^  and 
TTfc.  The  results  are  shown  in  Table  4.2. 

Table  4.2  Concentric  Tube  Heat  Exchanger  Outlet  Conditions  with  Typical  In- 


Flight  Inlet  Conditions 


Altitude 

Mach  Number 

Tta 

Pt4 

TTfc 

0  ft 

0 

1150  K 

1.09  MPa 

94.5% 

15,000  ft 

0.8 

1152  K 

0.93  MPa 

95.8% 

30,000  ft 

0.8 

1201  K 

0.49  MPa 

92.3% 

45,000  ft 

0.8 

1283  K 

0.25  MPa 

96.1% 

4.1.2  Trapezoidal  Radial  Fin  Geometry.  The  second  major  geometry  eval¬ 
uated  was  radially  oriented,  trapezoidal  fins.  The  intent  of  this  design  was  to  keep 
the  spacing  between  the  fins  constant  in  the  radial  direction  on  an  annular  heat 
exchanger.  For  this  purpose,  the  inner  radius  of  the  heat  exchanger  was  taken  to  be 
30.48  cm  (12  inches)  and  the  outer  radius  was  taken  to  be  50.8  cm  (20  inches).  With 
385  fins,  this  leaves  4.5  mm  spacing  throughout  and  fins  0.5  mm  thick  at  the  inner 
radius  and  3.8  mm  thick  at  the  outer  radius  (see  Figures  3.6  and  3.7  for  illustrations). 
Relatively  constant  heat  transfer  is  expected  in  the  radial  direction,  along  with  sig¬ 
nificantly  higher  pressure  loss  at  the  outer  radius,  just  due  to  pressure  drag  on  the 
thicker  fins.  The  total  temperature  at  the  outlet  was  expected  to  drop  with  increas¬ 
ing  fin  thickness,  due  to  the  higher  mass  fiow  through  the  channel  (since  the  density 
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and  velocity  are  held  constant  across  the  radial  direction,  the  thicker  fins  will  entail 
higher  mass  fiow,  as  they  capture  a  larger  area).  The  flow  in  this  three-dimensional 
area  was  approximated  with  two-dimensional  cuts  in  the  transverse  direction  (see 
Figures  3.7  and  3.8,  size  data  found  in  Table  3.1).  Data  from  these  cuts  was  linearly 
interpolated  to  form  a  distribution  in  the  radial  direction,  and  an  area  weighted  av¬ 
erage  total  temperature  and  pressure  were  calculated.  The  distributions  are  shown 
in  Figure  4.6  for  an  inlet  velocity  of  41  m/s  and  in  Figure  4.7  for  an  inlet  velocity 
of  29  m/s.  At  41  m/s  inlet  velocity,  the  average  Tta  'is  1106  K  (1990‘’R),  or  1040 
K  (1872‘’R)  corrected,  and  the  average  pressure  ratio  (TTft)  is  96.8%,  while  at  29  m/s 
inlet  velocity,  the  average  Tta  is  1212  K  (2181°R),  1130  K  (2033‘’R)  corrected,  and 
the  average  tt^  is  97.4%.  It  should  be  noted  that  the  total  pressure  and  temperature 
here  do  not  include  the  diffusion  or  compression  pressure  losses  of  the  inlet  or  outlet 
manifolds  and  possible  temperature  losses  due  to  cooling  air  injection.  From  the 
radial  inlet  and  outlet,  the  combined  estimated  pressure  ratio  of  the  inlet  and  outlet 
is  99.7%.  This  reduces  the  total  pressure  ratio  to  approximately  96.5%  for  41  m/s 
inlet  velocity  and  97.1%  for  29  m/s  inlet  velocity,  the  numbers  that  were  used  for 
thrust  calculations.  Again,  for  comparison  of  average  total  temperature  and  total 
pressure  ratio  with  other  geometries,  see  Table  4.6. 

As  with  the  concentric  tube  geometry,  heat  transfer  coefficients  were  evaluated 
to  ensure  that  they  stayed  within  reasonable  bounds.  The  variation  of  heat  flux  and 
heat  transfer  coefficient  along  the  center  of  the  trapezoidal  fin  is  shown  in  Figure 
4.8.  As  shown  in  Figure  4.8,  the  heat  transfer  coefficient  varied  from  about  1,000 

at  the  leading  edge  of  the  fin,  and  dropped  steadily  to  about  650  at  the 
trailing  edge  of  the  fin.  These  were  well  below  the  prediction  of  2,000  derived 
from  Kays  and  Crawford  (9:327).  Again,  the  data  matched  closely  with  predictions 
found  in  Hill  and  Peterson  (20:129),  and  since  the  other  cases  behaved  similarly,  the 
ANSYS®  (11)  data  was  taken  as  a  reasonable  approximation  of  reality. 
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Figure  4.6 


Figure  4.7 


|— Total  Temperature  — Total  Pressure] 


Radial  Variation  of  Total  Temperature  and  Total  Pressure  Across  Exit 
Plane,  Trapezoidal  Radial  Fin  Heat  Exchanger  at  41  m/s  Inlet  Velocity 
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Radial  Variation  of  Total  Temperature  and  Total  Pressure  Across  Exit 
Plane,  Trapezoidal  Radial  Fin  Heat  Exchanger  at  29  m/s  Inlet  Velocity 
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|— Heat  Flux - Film  CoeffTI 


Figure  4.8  Axial  Distribution  of  Heat  Flux  and  Heat  Transfer  Coefficient  at  Radial 
Center  of  Trapezoidal  Radial  Fins  at  41  m/s  Inlet  Velocity 

4.1.3  Rectangular  Radial  Fin  Geometry.  A  second  option  for  radial  fins 
was  the  rectangular  fin  with  variable  spacing.  This  used  two  annuli  of  2  mm  thick 
fins.  The  inner  annulus  had  and  inner  radius  of  30.48  cm  and  an  outer  radius  of 
40.64  cm.  With  344  fins,  this  gives  fin  spacing  from  3.57  mm  to  5.38  mm.  The 
outer  annulus  had  an  inner  radius  of  40.64  cm  and  an  outer  radius  of  50.8  cm. 
This  allowed  442  fins  with  spacing  from  3.78  mm  to  5.22  mm  (see  Figures  3.8, 
3.9,  and  3.10  for  illustrations  of  the  geometries  involved  in  this  case,  dimensions 
can  be  found  in  Table  3.2).  Each  annulus  was  approximated  with  three  transverse 
cuts  similar  to  those  used  for  the  trapezoidal  fins  (see  Figures  3.8  and  3.10,  and 
Table  3.2)  at  the  inner  radius,  outer  radius,  and  at  4.5  mm  spacing  (the  middle). 
The  data  was  linearly  interpolated  between  these  cuts  and  an  area-weighted  average 
was  calculated  from  the  interpolation.  The  distribution  of  total  pressure  and  total 
temperature  is  shown  in  Figure  4.9  for  41  m/s  inlet  velocity,  and  in  Figure  4.10  for 
29  m/s.  In  both  of  these  distributions  there  was  a  noticeable  discontinuity  in  total 
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temperature  and  total  pressure.  This  was  caused  by  the  transition  from  the  outer 
edge  of  the  inner  annulus  to  the  inner  edge  of  the  outer  annulus  (see  Figure  3.9  for 
a  schematic  view).  At  this  point  the  spacing  between  fins  changed  from  5.5  mm  to 
3.8  mm,  and  the  characteristics  of  the  flow  change  accordingly  (total  temperature 
rises  and  total  pressure  drops).  The  average  Tt4  41  m/s  inlet  velocity  was  1148 
K  (2066°R),  or  1076  K  (1936‘’R)  corrected,  and  the  average  pressure  ratio  (TTi)  at  41 
m/s  inlet  velocity  was  97.6  %.  At  29  m/s  inlet  velocity,  the  average  Tta  was  1221  K 
(2197‘’R),  1137  K  (2047°R)  corrected,  and  the  average  n  was  97.8%.  Again,  these 
total  pressures  and  temperatures  do  not  reflect  losses  in  the  inlet  or  outlet  manifolds, 
which  would  be  similar  to  the  estimates  for  the  trapezoidal  fins,  at  0.3%  total  pressure 
loss  (to  97.3%  for  the  41  m/s  inlet  velocity,  and  97.5%  for  the  29  m/s  inlet),  and 
an  indeterminate  total  temperature  loss  due  to  cooling  air.  For  comparison  of  these 
performance  parameters  to  those  for  other  geometries,  see  Table  4.6. 


"Total  Temperature  — Total  Pressure 


Figure  4.9  Radial  Variation  of  Total  Temperature  and  Total  Pressure  Across  Exit 
Plane,  Rectangular  Radial  Fin  Heat  Exchanger  at  41  m/s  Inlet  Velocity 
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Figure  4.10  Radial  Variation  of  Total  Temperature  and  Total  Pressure  Across  Exit 
Plane,  Rectangular  Radial  Fin  Heat  Exchanger  at  29  m/s  Inlet  Veloc¬ 
ity 


^.2  Heat  Exchanger  Performance  Trends 

By  recasting  some  of  the  data  taken  for  the  various  geometries,  some  basic 
trends  could  be  derived  as  to  which  parameters  (length,  spacing,  fin  thickness,  etc) 
affected  heat  exchanger  performance,  and  how  significantly.  By  rearranging  the  rect¬ 
angular  fin  data,  one  can  show  a  general  trend  in  the  performance  variation  due  to 
fin  spacing.  As  Figure  4.11  shows,  increasing  fin  spacing  decreased  the  total  tem¬ 
perature  available  and  increased  the  total  pressure  available.  However,  this  was  a 
very  small  change  in  total  pressure,  and  a  fairly  wide  variation  in  total  tempera¬ 
ture.  For  comparative  purposes,  a  single  run  was  made  with  2.0  mm  fins,  4.5  mm 
spacing,  and  508  mm  long  fins  rather  than  254  mm  fins,  while  keeping  all  bound¬ 
ary  conditions  identical  (at  41  m/s  inlet  velocity).  The  temperature  change  across 
the  heat  exchanger  rose  from  +433  K  (with  the  shorter  fins)  to  +  712  K  (with  the 
longer  fins)(+779°R  to  +1282°R),  while  the  pressure  loss  rose  from  1.8%  to  3.8%. 
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By  r6ca.sting  these  changes  into  the  form  of  a  Lj Dh  ratio,  and  taking  the  254  mm 
length,  9.0  mm  Dh  as  the  baseline,  it  can  be  shown  that  a  29%  increase  in  Lj Dh 
leads  to  a  42%  increase  in  temperature  change,  and  a  100%  increase  in  Lj  Dh  leads  to 
a  66%  increase  in  temperature  change.  While  this  shows  that  changes  to  fin  spacing 
or  fin  length  would  significantly  affect  the  outlet  temperatures,  the  fins  cannot  be 
simply  made  very  close  and  very  long  to  get  extremely  high  outlet  temperatures.  If 
the  solidity  (solidity  =  of  the  heat  exchanger  is  too  high,  the  mass  flow 

through  each  channel  will  increase  dramatically  {V  ^  partially  offsetting 

the  temperature  gains  (for  the  geometries  used  here  the  solidity  varies  from  30.8% 
for  the  dual  annulus  of  rectangular  fins,  to  32.4%  for  the  trapezoidal  fins,  with  the 
concentric  tubes  in  between  at  31.8%).  In  addition,  Fanno  Line  flow  (isothermal  flow 
in  a  constant  area  duct  with  friction)  and  Rayleigh  Line  flow  (frictionless  flow  in  a 
constant  area,  heated  duct)  show  that  there  is  a  limit  to  the  Lj  Dh  ratio  that  could 
be  used.  At  some  point,  from  the  combination  of  friction  and  heat  addition,  the 
Mach  number  will  reach  1,  at  which  point  the  flow  can  no  longer  accelerate,  and  the 
mass  flow  must  be  reduced.  The  increased  temperatures  will  tend  to  delay  this,  but 
will  lead  to  decreasing  total  pressure  also.  This  situation  would  reduce  the  efficiency 
of  the  engine  dramatically.  In  the  short  heat  exchangers  used  in  this  research,  this 
situation  is  not  present,  but  if  very  high  L/Dh  ratios  were  used,  this  could  be  a 
problem. 

The  data  for  both  the  rectangular  fins  and  the  trapezoidal  fins  demonstrate 
the  effect  of  lowering  inlet  velocity.  The  lower  velocity  impinging  on  the  fins  leads  to 
significantly  increased  outlet  temperatures.  Again,  this  is  due  to  the  mass  flow  drop¬ 
ping  faster  than  the  heat  transfer  coefficients.  A  further  demonstration  of  this  effect 
was  seen  in  the  altitude  variation  data  for  the  full-scale  concentric  heat  exchanger. 
This  data  was  recast  into  a  variation  in  Reynolds  number  based  on  hydraulic  di¬ 
ameter  {Reoh  =  ^^)-  For  this  comparison  the  velocity  is  derived  from  the  mass 
flow  into  the  burner  section  of  the  J-57  at  various  flight  conditions,  expanded  (in 
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Figure  4.11  Effects  of  Fin  Spacing  on  Corrected  Total  Temperature  and  Pressure 
at  Heat  Exchanger  Outlet  Plane,  at  29  m/s  Inlet  Velocity 

a  one-dimensional,  incompressible  calculation)  from  the  compressor  exit  to  an  area 
4  times  greater  (this  should  be,  on  the  average,  the  velocity  just  before  the  fins  or 
tubes  begin),  the  density  and  viscosity  are  derived  from  the  total  pressure  and  tem¬ 
perature  exiting  the  compressor,  and  the  hydraulic  diameter  of  the  tubes  is  twice 
the  spacing,  of  0.009  m.  The  effects  of  Reynolds  number  on  Tt4  and  are  shown 
in  Table  4.3  and  Figure  4.12. 


Table  4.3  Reynolds  Number  Effects  on  Tta  and  tt;, 


Reynolds  Number 

Tt4 

ATt 

^Pt 

TTb 

5.26  X  lO'^ 

1150  K 

-1-528  K 

-0.06  MPa 

94.5% 

4.08  X  10^ 

1152  K 

+b27K 

-0.04  MPa 

95.8% 

2.49  X  lO'^ 

1201  K 

+635  K 

-0.04  MPa 

92.2% 

1.29  X  10^ 

1283  K 

+747  K 

-0.01  MPa 

96.1  % 

4.2.1  Heat  Exchanger  Mass.  One  important  parameter  for  these  heat 
exchangers  will  be  the  mass  of  hafnium  contained  in  them.  This  will  drive  the 


4-14 


1300 


Reynolds  Number  (x10'^-4) 


Figure  4.12  Effects  of  Reynolds  Number  on  Available 

endurance  available,  as  well  as  the  structural  enhancements  required  to  support  the 
added  mass  of  the  heat  exchanger  in  the  engine.  One  area  of  some  uncertainty 
is  the  exact  density  of  the  hafnium  or  hafnium  alloy  as  it  would  be  used  in  the 
heat  exchanger.  Therefore,  a  mass  has  been  calculated  for  each  of  the  full-scale  heat 
exchangers  at  the  likely  extremes  of  the  density  range:  13,100  ^  and  9,000  This 
data  can  be  found  in  Table  4.4,  and  the  volume  of  the  fins  or  tubes  is  also  included. 
Note  that  the  trapezoidal  fin  heat  exchanger  mass  is  based  on  solid  hafnium  fins, 
although,  as  described  in  Section  3.3.10,  solid  isomer  fins  are  not  very  realistic  for 
this  geometry. 

Table  4.4  Mass  and  Solid  Volume  of  Heat  Exchangers 


Configuration 

Solid  Volume 

Mass  at 

p  =  13,100  ^ 

Mass  at 

P  =  9,000  ^ 

Concentric  Tubes 
Trapezoidal  Fins 
Rectangular  Fins 

0.042  rri^ 
0.043 

0.041 

549  kg 

560  kg 

531  kg 

377  kg 

386  kg 

365  kg 

4-15 


Even  on  an  engine  as  large  as  the  J-57  (1,757  kg  or  3,870  lb),  this  is  a  signif¬ 
icant  addition  of  mass  (365-560  kg).  To  this  we  must  also  add  several  hundred  or 
thousand  pounds  of  radiation  shielding  either  at  the  source  (the  engine)  or  around 
sensitive  components,  or  both.  This  weight  gain  is  offset,  though,  by  the  loss  of 
several  thousand  pounds  of  fuel  weight.  These  factors  would  be  driven  by  the  ap¬ 
plication,  but  it  is  obvious  that  the  mass  increase  of  the  engine  will  require  some 
significant  structural  changes.  While  the  exact  performance  effects  of  this  additional 
mass  cannot  be  assessed  without  information  about  the  application,  it  will  not  be 
negligible. 

4.. 3  Engine  Performance  Estimates 

4.3.1  Static  Sea- Level  Performance.  Performance  estimates  for  the  J-57 
turbojet  with  these  heat  exchangers  replacing  the  conventional  combustion  chamber 
were  accomplished  using  the  ONX  software  provided  with  Mattingly,  Heiser,  and 
Daley’s  book  Aircraft  Engine  Design.  This  software  enables  the  user  to  change 
individual  engine  parameters  to  get  an  on-design  thrust  estimate  for  the  engine 
(4:97-133;23). 

To  evaluate  the  performance  of  the  J-57  with  the  heat  exchangers  in  place, 
it  was  first  necessary  to  calculate  the  other  parameters  for  the  engine,  such  as: 
compressor  pressure  ratio,  compressor  and  turbine  efficiencies,  and  so  forth.  These 
are  shown  in  Table  4.5.  By  taking  these  parameters  to  be  fixed,  the  thrust  of  the 
engine  can  be  estimated  with  the  heat  exchangers  in  place  by  changing  Tn  and  tti,. 
The  results  of  these  calculations  are  found  in  Table  4.6.  Please  note  that  the  pressure 
ratios  for  the  radial  fin  geometries  include  estimated  pressure  losses  for  the  inlet  and 
outlet  manifolds,  derived  from  the  radial  inlet  and  outlet  models. 

All  three  of  the  full-size  geometries  exceeded  the  thrust  value  of  the  conven¬ 
tional  J-57  turbojet  (10,344  IbT,  as  calculated  using  ONX) (23).  The  reader  should 
note  that  the  higher  velocity  (41  m/s)  runs  with  the  radial  fin  geometries  would 
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Table  4.5  Critical  Engine  Parameters  for  the  J-57  Turbojet 


Calculated  Values 

7rc/(Low  Pressure  Compressor  Pressure  Ratio) 

3.673 

7rc/i(High  Pressure  Compressor  Pressure  Ratio) 

3.092 

77c/ (Low  Pressure  Compressor  Poly  tropic  Efficiency) 

0.885 

77c/i(High  Pressure  Compressor  Polytropic  Efficiency) 

0.923 

77t  (Turbine  Polytropic  Efficiency) 

0.938 

Assumed  Values 

7c(Cold  Ratio  of  Specific  Heats) 

1.4 

7/i(Hot  Ratio  of  Specific  Heats) 

1.3 

(Maximum  Diffuser  Pressure  Ratio) 

0.97 

T^nmax  (Maximum  Nozzle  Pressure  Ratio) 

0.98 

77n,/ (Mechanical  Efficiency,  Low  Pressure  Spool) 

0.98 

77m/i  (Mechanical  Efficiency,  High  Pressure  Spool) 

0.98 

Bleed  Air  Percentage 

1.0% 

Cooling  Air  Take-off  #1 

5.0% 

Cooling  Air  Take-off  #2 

5.0% 

Power  Take-off 

250  kW 

?fi(Engine  mass  air  flow) 

167  Ibm/s 

Ta 


ble  4.6  Estimated  J-57  Sea  Level  Static  Thrust  with  Heat  Exchangers  Instal 


Heat  Exchanger  Type 

Tta 

(K) 

Thrust 

(IbT) 

%  of  J-57 

0.5  mm  Concentric  Tubes 

1027  K 

96.8% 

8,964 

86.6% 

2.0  mm  Concentric  Tubes 

986  K 

97.0% 

8,301 

80.2% 

0.5  mm  Staggered  Concentric  Tubes 

1042  K 

95.0% 

9,105 

88.0% 

Full-Scale  Concentric  Tubes 

1150  K 

94.5% 

10,619 

102.7% 

Trapezoidal  Radial  Fins  (41  m/s) 

1040  K 

96.5% 

9,147 

88.4% 

Trapezoidal  Radial  Fins  (29  m/s) 

1130  K 

97.1% 

10,488 

101.4% 

Rectangular  Radial  Fins  (41  m/s) 

1076  K 

97.3% 

9,736 

94.1% 

Rectangular  Radial  Fins  (29  m/s) 

1137  K 

97.5% 

10,607 

102.5% 

led 
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represent  a  higher  mass  flow  rate  through  the  heat  exchanger.  An  increased  mass 
flow  rate  would  alter  the  performance  from  that  listed  here,  but  would  not  represent 
the  effect  of  installing  these  heat  exchangers  in  an  off-the-shelf  engine.  The  higher 
velocity  data  would  also  represent  the  effect  of  installing  a  smaller  (radially)  heat 
exchanger  to  accommodate  space  restrictions.  It  should  also  be  noted  that  these 
static,  sea  level  thrust  values  only  represent  part  of  the  overall  performance  picture. 
They  do  not  represent  the  variation  in  heat  exchanger,  and  engine,  performance 
which  would  occur  due  to  altitude  and  ofif-design  operation  conditions. 

4.S.2  Performance  at  Other  Flight  Conditions.  From  the  data  for  the 
full-scale  concentric  tube  geometry,  the  flight  performance  of  the  heat  exchanger 
can  be  seen.  The  heat  exchanger  proved  capable  of  heating  the  flow  to  higher 
temperatures  at  higher  altitudes.  The  J-57,  though,  cannot  use  this  excess  capacity. 
As  the  aircraft  climbs,  the  engine  needs  less  total  temperature  capability  to  allow 
the  compressor  to  operate  at  its  limiting  pressure  ratio.  Therefore,  as  the  aircraft 
climbs,  the  requirement  for  turbine  inlet  temperature  drops.  This  trend  is  shown 
in  Figure  4.13.  Since  the  heat  exchanger  proved  able  to  produce  more  turbine  inlet 
temperature  than  needed,  the  requirement  could  be  matched  by  lowering  the  average 
wall  temperature  until  the  7r4  requirement  was  met.  Under  these  conditions,  the 
engine  with  heat  exchangers  installed  should  very  nearly  match  the  baseline  J-57  s 
performance  under  flight  conditions. 

4.4  Grid  Independence 

Due  to  the  size  of  the  concentric  models,  and  the  fact  that  spacing  was  nearly 
identical  in  the  radial  fin  models,  grid  independence  was  checked  with  a  radial  fin 
model,  with  0.5  mm  thick  fins  and  4.5  mm  spacing.  The  results  showed  a  1  degree 
increase  in  outlet  temperature  (about  a  0.2%  change)  and  a  0.2%  increase  in  outlet 
total  pressure.  Therefore,  the  data  was  largely  grid  independent. 
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V.  Conclusions  and  Recommendations 

This  chapter  contains  conclusions  based  on  the  current  work  and  recommen¬ 
dations  concerning  directions  for  future  research  in  this  area. 

In  this  research,  three  heat  exchanger  configurations  (concentric  tubes,  radial 
fins  with  constant  spacing,  and  radial  fins  with  constant  thickness)  were  evaluated  for 
their  performance  in  a  gas-turbine  environment.  The  performance  ratings  used  were 
the  exit  total  temperature  and  the  total  pressure  loss.  In  a  gas  turbine  engine  these 
translate  to  the  turbine  inlet  temperature,  Tt^,  and  burner  pressure  ratio,  tt^.  The 
intended  application  of  this  research  would  be  the  use  of  a  radioactive  isomer  whose 
rapid  decay  can  be  triggered  and  controlled,  to  power  a  high  endurance  jet  aircraft. 
The  research  was  based  on  the  following  assumptions:  1)  That  the  decay  process  can 
be  controlled.  2)  That  the  isomer,  and  whatever  materials  it  may  be  alloyed  with  or 
contained  in,  have  similar  material  properties  to  pure  metallic  hafnium  (especially 
the  melting  point).  3)  That  the  heated  surfaces  of  the  heat  exchanger  would  be 
manufactured  from  the  isomer  and  its  associated  materials.  Assumptions  in  these 
areas  were  necessary  to  begin  analysis,  and  these  were  judged  reasonable. 

Heat  exchanger  geometries  were  intended  to  fit  within  the  envelope  of  the 
combustion  chamber  section  of  a  J-57  turbojet  engine.  The  J-57  was  also  used 
as  a  source  for  entrance  conditions  to  the  heat  exchangers,  and  as  a  baseline  for 
performance  comparisons. 

5.1  Conclusions 

The  following  conclusions  are  presented  based  on  the  data  in  this  research,  and 
the  entering  assumptions. 

1.  Any  of  these  three  configurations,  if  it  can  be  manufactured,  and  the  heat 
generation  rates  can  be  controlled  at  a  close  to  uniform  level,  should  be  a  feasible 
replacement  for  the  combustion  chamber  in  a  turbojet.  This  feasibility  is  based  only 
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on  the  thermodynamics  and  fluid  mechanics  of  the  problem,  and  does  not  consider 
application  driven  or  operational  considerations.  The  manufacturability  of  these 
heat  exchangers  is  largely  a  machining  question,  and  a  question  of  the  availability  of 
large  quantities  of  the  isomer. 

2.  Engine  operations  at  higher  altitudes  should  be  easier  than  at  sea  level.  In 
general,  the  temperature  requirements  will  be  lower,  while  the  capability  of  the  heat 
exchanger  to  raise  the  temperature  of  the  flow  is  increasing.  In  addition,  the  lower 
mass  flow  rates  at  high  altitudes  will  reduce  the  energy  and  radiation  output  of  the 
heat  exchanger,  extending  endurance  and  limiting  the  radiation  exposure  of  sensitive 
components  of  the  aircraft. 

3.  The  extraordinary  endurance  oflFered  by  this  heat  source  could  drastically 
change  the  operating  paradigms  for  many  missions.  Overall,  this  research  showed 
that  an  isomer  heat  source  may  be  feasibly  harnessed  to  power  a  turbojet  (or  tur¬ 
bofan)  engine,  if  the  radiation  problem  can  be  dealt  with. 

5.2  Recommendations 

There  are  several  directions  that  could  be  taken  from  this  work.  Some  of  these 
are  listed  below,  in  no  particular  order  of  precedence. 

1.  Optimization  of  the  heat  exchanger  geometry  for  a  given  airframe  and 
engine  geometry.  For  example,  given  hard  limits  on  the  turbine  inlet  temperature 
and  pressure  capacity  of  an  engine,  design  a  heat  exchanger  geometry  for  maximum 
thrust,  or  maximum  endurance.  Obviously,  this  would  require  some  knowledge  of 
the  proposed  mission  for  the  aircraft.  An  engine  design  could  even  be  optimized  for 
minimal  radiation,  which  would  probably  lead  to  engine  designs  similar  to  those  for 
maximum  fuel  efficiency,  with  low  hot  mass  ffows,  and  high  bypass  ratios.  To  take 
this  path  would  require  some  specification  of  airframe  and  engine  data. 

2.  Some  consideration  must  be  give  to  aeroelasticity  in  this  heat  exchanger. 
The  minimum  stiffness  of  the  fins  (calculated  from  the  material  properties  and  ge- 
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ometry)  which  can  be  used  should  be  determined.  From  the  material  properties  for 
hafnium  listed  in  the  ASM  Handbook  (27:665),  this  could  be  quite  thin,  but  this 
data  did  not  include  temperature  variations  of  these  properties.  Closely  related  is 
the  resistance  to  corrosion  or  erosion  of  the  hafnium  in  this  high-temperature  envi¬ 
ronment  with  constant  airflow,  probably  carrying  particulates  of  some  type.  Detailed 
studies  of  the  material  behavior  at  high  temperatures,  up  to  the  melting  point,  may 
be  required. 

3.  Structural  considerations  also  need  to  be  addressed.  Support  of  this  very 
high  temperature,  very  heavy  heat  exchanger  should  be  addressed.  The  use  of  some 
type  of  ceramic  supports,  separators,  or  insulators  could  be  evalutated.  Air  cooling 
requirements  for  the  structural  members  of  the  heat  exchanger  could  be  included. 
These  are  some  of  the  structural  issues  which  would  require  analysis  before  an  op¬ 
erational  engine  could  be  built. 

4.  The  use  of  still  thicker  fins  could  be  evaluated,  with  the  aim  of  using  hollow 
fins  made  of  some  highly  refractory  material,  such  as  titanium  carbide  or  hafnium 
carbide,  to  contain  liquid  hafnium  isomer.  This  might  allow  wall  temperatures  of 
3,000  K  or  even  higher.  This  could  allow  much  higher  heat  transfer  or  more  compact 
heat  exchangers.  Obviously,  the  chemical  reactivity  of  the  hafnium  at  that  point 
would  need  to  be  evaluated,  as  well  as  the  properties  of  the  fin  casing  material  at 
extremely  high  temperatures. 

5.  Given  that  the  use  of  this  isomer  in  a  gas-turbine  engine  is  feasible,  a  de¬ 
tailed  study  of  the  radiation  shielding  requirements  may  be  in  order.  This  would 
probably  include:  the  exact  spectrum  of  radiation  emitted  by  the  heat  exchanger 
in  full-throttle  operation,  some  determination  of  the  maximum  radiation  dose  per¬ 
missible  to  the  most  sensitive  components  (the  pilot  if  the  aircraft  is  manned),  and, 
given  an  airframe  geometry  (engine  location  relative  to  sensitive  components)  and 
some  material  limitations,  design  radiation  shielding  to  limit  radiation  exposure  to 
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acceptable  level.  This  area  of  study  is  probably  best  left  to  physicists  or  nuclear 
engineers,  who  have  a  background  in  radiation  shielding. 

6.  Some  study  could  be  made  of  the  operational  considerations  this  power 
source  would  entail.  This  could  include,  but  is  not  limited  to.  considerations  for 
maintenance  personnel  who  have  to  work  on  the  engine  or  aircraft;  operational  safety 
considerations  including  take-off  and  landing  radiation  hazards  to  personnel  at  the 
airfield;  disposal  of  used  heat  exchanger  fins;  crash  safety  of  the  heat  exchangers; 
and  political  considerations  of  operating  this  type  of  aircraft  over  populated  areas 
or  in  foreign  airspace.  All  of  these  considerations  would  need  to  be  addressed  before 
this  type  of  system  could  be  operational. 


Appendix  A.  Temperature  Corrections  for  Variable  Specific  Heats 

Since  ANSYS®  (11)  cannot  use  a  variable  specific  heat  with  a  compressible 
flow,  it  was  necessary  to  correct  the  temperature  data  at  the  end  of  the  process. 

This  correction  was  based  on  two  premises:  first,  that  the  heat  transfer  data  is 
not  significantly  affected  by  the  temperature  errors  introduced  by  the  fixed  specific 
heat;  and  second;  that  the  variation  in  specific  heats  can  be  adequately  represented 
by  the  algebraic  average  of  the  specific  heats  at  the  high  and  low  temperatures 
in  the  process.  As  for  the  first  assumption,  the  heat  transfer  coefficients  are  largely 
determined  by  the  Reynolds  number  of  the  flow  at  the  specific  location.  Temperature 
has  direct  effects  on  two  parts  of  the  Reynolds  number:  the  density  and  the  viscosity. 
The  density,  however,  is  part  of  a  product  with  velocity.  Since  the  mass  flow  is 
constant,  therefore,  the  Reynolds  number  is  unchanged  by  fluctuations  in  density. 
The  temperature  changes  in  viscosity  are  not  similarly  negated,  but  are  fairly  small 
in  the  range  of  the  temperature  change  at  any  given  location  in  the  flow  due  to  the 
specific  heat  variation  (less  than  100  K).  For  the  second  assumption,  a  check  of  the 
variation  of  flow  enthalpy  with  temperature  shows  that  the  average  specific  heats 
between  two  points  are  determined  by: 


Cp  — 


ATi,2 
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This  is  very  close  (within  <1%)  to  the  algebraic  average  of  the  Cp  for  the  end 
points.  Therefore  the  outlet  temperatures  were  corrected  using  the  following  method. 

Step  1:  The  enthalpy  change  is  computed  from: 


Ah  =  CpANSYsi^Tout  -  TXin)  (■^•2) 
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Step  2:  The  Cp  for  the  inlet  temperature  is  computed  from  a  correlation  found 
in  Wilson  and  Korakianitis  (32:575).  This  correlation  is  of  the  form: 
with  Ci  given  for  ranges  of  temperatures. 

Step  3:  The  Cp  for  the  initial  outlet  temperature  is  computed  from  the  appro¬ 
priate  correlation  in  Wilson  and  Korakianitis,  and  an  average  Cp  is  calculated. 


ow  ^high 
2 


(A.3) 


Step  4:  The  new  outlet  temperature  is  generated  from  the  following  equation: 


Ah  ^ 

=  ;; —  +  Ttz 

^avg 


(A.4) 


Step  5:  A  new  Cp^^^  is  computed,  based  on  the  new  Tn,  and  the  process  repeats 
until  a  converged  value  is  established. 
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Appendix  B.  Building  a  Model  in  ANSYl^ 

This  Appendix  covers  the  basics  of  building  a  simple  heat  exchanger  model  in 
ANSYS®  5.6.1(Faculty  Research  Release)  (11).  These  instructions  were  based  on 
the  methods,  software,  and  hardware  used  in  this  research.  ANSYS®  5.6.1  was  run 
on  a  700  MHz  Pentium  III®  PC,  with  384  Mb  RAM,  a  19”  monitor,  and  2.0  Gb 
hard  drive  (if  purchasing  hardware  to  run  this  software,  maximize  RAM  and  storage, 
minimize  fluff  such  as  speakers  and  sound  cards.  A  CD  ROM  drive  is  mandatory  for 
installation,  large  capacity  backup  drives  would  be  useful.) 

Step  1.  Start  ANSYS®.  (in  Windows®:  PROGRAMS,  ANSYS®,  INTERAC¬ 
TIVE) 

Step  2.  Enter  Initial  Jobname  (ex:  heatexchangerl),  press  RUN. (If  you  specify 
an  existing  file,  it  will  open  that  file.  If  you  do  not  want  to  overwrite  the  old  file, 
you  should  press  RESUME-DB  as  soon  as  ANSYS®  starts.)  This  jobname  will  be 
the  name  of  all  files  associated  with  this  model. 

Step  3.  Go  to  PREFERENCES,  press  FLOTRAN®  CFD,  press  OK.  This  tells 
ANSYS  that  you  want  to  run  a  CFD  model.  Note  that  no  other  types  of  elements 
can  be  present  in  CFD  models. 

Step  4.  Go  to  PREPROCESSOR,  ELEMENT  TYPE,  ADD/EDIT/DELETE, 
ADD,  select  2D  FLOTRAN®  141,  OK,  CLOSE.  Close  ELEMENT  TYPE  Menu. 
This  tells  ANSYS  that  you  will  be  running  a  2-D  analysis.  Again,  only  FLOTRAN 
element  types  are  permitted  in  a  CFD  analysis. 

Step  5.  Go  to  CREATE  Submenu.  Select  AREA,  RECTANGLE,  BY  DIMEN¬ 
SIONS.  Enter  the  following  rectangles: 
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XI 

X2 

Y1 

Y2 

0 

0.0254 

0.0127 

0.0381 

0.0762 

0.3302 

0 

0.054 

0.381 

0.4064 

0.0127 

0.0381 

0.4064 

0.75 

0.0127 

0.0381 

This  creates  the  entrance,  body,  and  exit  of  the  heat  exchanger.  The  dimen¬ 
sions  are  all  in  meters. 

Step  6.  Select  LINE,  LINES,  TAN  TO  2  LINES,  connect  the  top  and  bottom 
of  the  first  rectangle  to  the  top  and  bottom  of  the  second,  and  the  top  and  bottom 
of  the  second  to  the  top  and  bottom  of  the  third,  following  the  directions  on  the 
screen. 

Step  7.  Select  AREA,  ARBITRARY,  BY  LINES.  Select  each  of  the  four  lines 
making  up  the  expanding  and  contracting  areas  on  either  side  of  the  second  rectangle. 
This  makes  these  into  ANSYS®  areas.  Steps  6  and  7  create  smooth  transitions  from 
the  inlet  to  the  body  and  the  body  to  the  outlet. 

Step  8.  Close  CREATE  menu.  Open  OPERATE  menu,  select  ADD,  AREAS. 
Select  all  areas.  This  makes  all  six  areas  into  one  large  area.  This  completes  the 
envelope  of  the  heat  exchanger. 

Step  9.  Close  OPERATE  menu.  Open  CREATE  menu.  Select  AREA,  RECT¬ 
ANGLE,  BY  DIMENSIONS.  Create  the  following  rectangles: 
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XI 

X2 

Y1 

Y2 

0.0762 

3302 

0.0045 

0.0055 

0.0762 

3302 

0.01 

0.011 

0.0762 

3302 

0.0155 

0.0165 

0.0762 

3302 

0.021 

0.022 

0.0762 

3302 

0.0265 

0.0275 

0.0762 

3302 

0.032 

0.033 

0.0762 

3302 

0.0375 

0.0385 

0.0762 

3302 

0.043 

0.044 

0.0762 

3302 

0.0485 

0.0495 

This  creates  the  nine  fins  or  tubes  of  this  heat  exchanger.  Each  is  254  mm  long  and 
1  mm  thick. 

Step  10.  Close  CREATE  menu.  Open  OPERATE  menu.  Select  SUBTRACT. 
Select  the  large  envelope  area  (this  is  the  area  we  subtract  from).  Enter.  Select  the 
nine  small  rectangles  just  created  (these  will  be  subtracted).  Enter.  Now  the  area  of 
the  heat  exchanger  has  nine  voids  where  the  fins  are,  and  flow  is  not  allowed.  The 
picture  on  your  screen  should  look  very  similar  to  the  compact  concentric  geometry 
shown  in  Chapter  3,  Figure  3.2.  Close  OPERATE  menu. 

Step  11.  Open  MESHTOOL.  Select  LINE-SET.  Set  all  lines  to  equal  spacing 
at  about  0.6  mm  per  division  (Number  of  Divisions  =  1.6  •  line  length  in  mm, 
all  dimensions  were  entered  in  m).  Select  Tri,  Free,  MESH,  Select  area  on  screen, 
ENTER.  ANSYS®  will  now  generate  a  grid/mesh.  Select  SAVE-DB  at  top  of  screen. 
The  grid  should  now  be  present,  and  will  probably  be  fine  enough  to  meet  y+ 
requirements,  which  we  will  check  later. 

Step  12.  Open  FLOTRAN  SET  UP.  Select  settings  as  follows. 

SOLUTION  OPTIONS  Select:  Transient 

Solve  Flow  Equations 
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Thermal 

Turbulent 

Compressible 

No  for  all  other  options 

EXECUTION  OPTIONS  Select  OK  (accept  Defaults  on  first  menu 

5000  Global  Iterations 
Initial  time  step  value  =  0 
Number  of  time  steps  =  5000 
Global  Iterations  per  time  step  =  20 
accept  all  other  defaults 

ADDITIONAL  OUT  Select:  RFL  Out  Derived 

YPLU 

TAUW 

FLUID  PROPERTIES  Select:  AIR-SI  for  all  properties 

FLOW  ENVIRONMENT  Select:  Ref  Conditions 

REF  CONDITIONS  Select:  Reference  Pressure  =  1151000  (Pa) 

Ratio  of  Cp/Cv  =  1.35 
Nominal  Temperature  =  622  K 
Total  Temperature  =  622  K 
Bulk  Temperature  =  622  K 
Temp  offset  from  absolute  zero  =  0 
CFD  SOLVER  CONTR  Select:  VX  Solver  CFD  =  TDMA 

VY  Solver  CFD  =  TDMA 
TEMP  Solver  CFD  =  TDMA 
PRES  Solver  CFD  =  PCRS 
ENKE  Solver  CFD  =  PGMR 
ENDS  Solver  CFD  =  PGMR 


TURBULENCE  Select:  Turbulence  Model 
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ADVECTION 


Shi-Zhu-Lumley 

Select:  Momentum  Equations  =  MSU 
Turbulent  Equations  =  SUPG 
Compressible  Pressure  Equations  =  MSU 
Energy  Equation  =  MSU 


Most  of  these  entries  should  be  self-explanatory.  These  will  set  up  ANSYS/Flotran 
to  perform  the  analysis  for  a  heat  exchanger. 

Step  13.  Open  LOADS,  APPLY.  Apply  Constant  Temperature  to  all  lines 
making  up  the  fins  (2400).  Apply  Constant  Heat  Flux  (0)  to  all  other  lines,  save 
inlet  and  outlet.  Apply  Velocity  (X-velocity  =  0,  Y-velocity  =  0)  to  all  lines  save 
inlet  and  outlet.  Apply  Velocity  (X-velocity  =  75,  Y- Velocity  =  0)  to  inlet.  Apply 
Pressure  DOF  (0)  to  outlet.  Apply  INITIAL  CONDITION,  SELECT  ALL,  VX  = 

75.  Close  Loads.  Close  PREPROCESSOR.  Select  SAVE-DB.  These  settings  will 
prohibit  flow  crossing  solid  boundaries  of  the  model,  impose  shear  forces  along  the 
solid  boundaries,  impose  a  heat  transfer  load  on  the  heated  isomer  fins  of  the  model, 
and  impose  an  insulated  wall  condition  on  the  unheated  walls  of  the  model.  The  0 
(relative)  pressure  at  the  outlet  is  required  for  subsonic  flow,  and  allows  calculation 
of  pressures  upstream  relative  to  the  outlet  condition. 

Step  14.  Open  SOLUTION.  Select  FLOTRAN®  SET  UP.  Select  FLOCHECK, 
INITIALIZE  ALL,  ENTER.  Select  SAVE-DB.  This  applies  the  initial  conditions  and 
boundary  conditions  to  all  of  the  elements  in  the  model. 

Step  15.  Select  RUN  FLOTRAN®. 

Step  16.  When  complete  (or  Program  may  be  stopped  using  DONE  key  on 
screen),  Select  SAVE-DB.  Close  SOLUTION.  Open  GENERAL  POSTPROCES- 
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SOR.  Select  READ  RESULTS,  LAST  SET.  This  enables  the  postprocessor  to  use 
the  current  solution  output  to  create  all  of  the  output  displays  later. 

Step  17.  Select  PLOT  RESULTS,  NODAL  SOLUTION.  This  will  plot  contours 
of  velocity,  temperature,  pressure,  and  so  forth.  Explore  these  at  will.  Return 
to  GENERAL  POSTPROCESSOR  menu.  These  figures  (or  any  other  image  in 
the  ANSYS®  Graphics  window)  can  be  saved  through  PLOTCTRLS,  CAPTURE 
IMAGE.  This  will  create  a  bitmap  (.BMP)  file  which  can  be  saved  in  the  filename 

of  choice. 

Step  18.  Select  PATH  OPERATIONS.  Select  DEFINE  PATH,  BY  NODES, 
Select  the  top  and  bottom  of  the  inlet,  use  20  divisions,  default  data  sets.  Return 
to  PATH  OPERATIONS.  Select  MAP  ONTO  PATH.  Select  VX,  TEMP,  PRES, 
TTOT,  PTOT.  Select  PLOT  PATH  ITEM.  This  will  generate  plots  of  the  velocity, 
pressure,  etc,  along  the  path.  Explore  at  will.  Select  LIST  PATH  ITEM.  Select  VX, 
TEMP,  PRES,  TTOT,  PTOT.  This  will  create  a  listing  of  the  variables  and  distance 
along  the  path  specified.  The  path  operations  can  be  repeated  for  any  path  desired. 
Save  Path  item  listings  (FILE.  SAVE  AS),  and  they  can  be  opened  using  EXCEL® 
or  some  other  tool  for  further  analysis. 

These  instructions  should  get  a  new  user  started.  There  are  many  subtleties 
to  using  ANSYS®,  and  only  experience  will  uncover  them  all.  As  with  any  other 
computer  tool,  SAVE  OFTEN.  ANSYS®  does  consume  a  significant  amount  of  hard- 
drive  space,  so  plan  accordingly. 

Notes: 

1.  If  element  counts  are  above  100,000,  runs  will  be  quite  slow,  and  this 
version  of  ANSYS®  cannot  handle  over  250,000  elements.  A  50,000  element  model 
will  occupy  about  100  Mb  of  disk  space  for  model,  solution,  and  associated  files. 
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2.  Transient  solutions  seem  more  stable  than  steady-state  solutions.  It  may 
be  possible  to  run  the  problem  transient  for  a  period  and  then  solve  as  steady  state 
(not  always). 

3.  Remember  that  compressible  flows  must  have  flxed  speciflc  heats.  Try  to 
use  a  good  guess  at  speciflc  heat,  and  then  correct  the  output  if  necessary. 

4.  Remember  that  all  pressure  outputs  are  relative  to  the  reference  pressure 
that  was  entered  in  (PREPROCESSOR,  FLOTRAN  SET  UP,  FLOW  ENVIRON¬ 
MENT).  When  using  these  temperatures  in  other  programs  you  will  need  to  add  the 
reference  pressure  back  in. 
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